Topological Defects and Gapless Modes in Insulators and Superconductors 



Jeffrey C.Y. Teo and C.L. Kane 
Dept. of Physics and Astronomy, University of Pennsylvania, Philadelphia, PA 19104 

We develop a unified framework to classify topological defects in insulators and superconduc- 
tors described by spatially modulated Bloch and Bogoliubov de Gennes Hamiltonians. We consider 
Hamiltonians 'H(k, r) that vary slowly with adiabatic parameters r surrounding the defect and belong 
to any of the ten symmetry classes defined by time reversal symmetry and particle-hole symmetry. 
The topological classes for such defects are identified, and explicit formulas for the topological invari- 
ants are presented. We introduce a generalization of the bulk-boundary correspondence that relates 
the topological classes to defect Hamiltonians to the presence of protected gapless modes at the de- 
fect. Many examples of line and point defects in three dimensional systems will be discussed. These 
can host one dimensional chiral Dirac fermions, helical Dirac fermions, chiral Majorana fermions 
and helical Majorana fermions, as well as zero dimensional chiral and Majorana zero modes. This 
approach can also be used to classify temporal pumping cycles, such as the Thouless charge pump, 
as well as a fermion parity pump, which is related to the Ising non-Abelian statistics of defects that 
support Majorana zero modes. 

PACS numbers: 73.20.-r,73.43.-f, 71.10.Pm, 74.45.-l-c 



I. INTRODUCTION 

The classification of electronic phases according to 
topological invariants is a powerful tool for understand- 
ing and predicting the behavior of matter. This approach 
was pioneered by Thouless, et alfi(TKNN), who identi- 
fied the integer topological invariant characterizing the 
two dimensional (2D) integer quantum Hall state. The 
TKNN invariant n gives the Hall conductivity a^y = 
ne^/h and characterizes the Bloch Hamiltonian 'H(k), 
defined as a function of k in the magnetic Brillouin zone. 
It may be expressed as the first Chern number associ- 
ated with the Bloch wavefunctions of the occupied states. 
A fundamental consequence of this topological classifica- 
tion is the hulk-boundary correspondence^ which relates 
the topological class of the bulk system to the number of 
gapless chiral fermion edge states on the sample bound- 
ary. 

Recent interest in topological states^"— has been stim- 
ulated by the realization that the combination of time 
reversal symmetry and the spin orbit interaction can 
lead to topological insulating electronic phases^"— and by 
the predictionii"— and observationi^"— of these phases 
in real materials. A topological insulator is a two or 
three dimensional material with a bulk energy gap that 
has gapless modes on the edge or surface that are pro- 
tected by time reversal symmetry. The bulk boundary 
correspondence relates these modes to a Z2 topologi- 
cal invariant characterizing time reversal invariant Bloch 
Hamiltonians. Signatures of these protected boundary 
modes have been observed in transport experiments on 
2D HgCdTe quantum wellsi^"— and in photoemission 
and STM experiments on 3D crystals of Bii_a;Sbar^"— , 
BiaSea^, Bhle^^^^^^ and SbaTe^^. Topological insu- 
lator behavior has also been predicted in other classes of 
materials with strong spin orbit interactions^i"— . 

Superconductors, described within a Bogoliubov de 
Gennes (BdG) framework can similarly be classi- 



fied topologicallji^— . The Bloch-BdG Hamiltonian 
^sdG(k) has a structure similar to an ordinary Bloch 
Hamiltonian, except that it has an exact particle-hole 
symmetry that reflects the particle-hole redundancy in- 
herent to the BdG theory. Topological superconductors 
are also characterized by gapless boundary modes. How- 
ever, due to the particle-hole redundancy, the boundary 
excitations are Majorana fermions. The simplest model 
topological superconductor is a weakly paired spinless 
p wave superconductor in ID^, which has zero energy 
Majorana bound states at its ends. In 2D, a weakly 
paired Px + ipy superconductor has a chiral Majorana 
edge stated. Sr2Ru04 is believed to exhibit a triplet 
Px + iPy stated. The spin degeneracy, however, leads to 
a doubling of the Majorana edge states. Though undou- 
blcd topological superconductors remain to be discovered 
experimentally, superfiuid '^He B is a related topological 
phase^i^'^i^i^ and is predicted to exhibit 2D gapless 
Majorana modes on its surface. Related ideas have also 
been used to topologically classify Fermi surfaces^. 

Topological insulators and superconductors fit to- 
gether into an elegant mathematical framework that gen- 
eralizes the above classifications'^*^. The topological 
classification of a general Bloch or BdG theory is specified 
by the dimension d and the 10 Altland Zirnbauer sym- 
metry classes^ characterizing the presence or absence of 
particle-hole, time reversal and/or chiral symmetry. The 
topological classifications, given by Z, Z2 or show a 
regular pattern as a function of symmetry class and d, 
and can be arranged into a periodic table of topological 
insulators and superconductors. Each non trivial entry 
in the table is predicted, via the bulk-boundary corre- 
spondence, to have gapless boundary states. 

Topologically protected zero modes and gapless states 
can also occur at topological defects, and have deep im- 
plications in both field theory and condensed matter 
physics^i^"— . A simple example is the zero energy 
Majorana mode that occurs at a vortex in a, px -\- ipy 
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FIG. 1: Topological defects characterized by a _D parameter 
family of d dimensional Bloch-BdG Hamiltonians. Line de- 
fects correspond to d — _D = 2, while point defects correspond 
to d — D = 1. Temporal cycles for point defects correspond 
to d - D = 0. 



superconductor—. Similar Majorana bound states can 
be engineered using three dimensional heterostructures 
that combine ordinary superconductors and topological 
insulators^, as well as semiconductor structures that 
combine superconductivity, magnetism and strong spin 
orbit interactions^"—. Recently, we showed that the ex- 
istence of a Majorana bound state at a point defect in a 
three dimensional Bogoliubov de Gennes theory is related 
to a Z2 topological invariant that characterizes a family 
of Bogoliubov de Gennes Hamiltonians H. BdG (k, r) de- 
fined for r on a surface surrounding the defecli^. This 
suggests that a more general formulation of topological 
defects and their corresponding gapless modes should be 
possible. 

In this paper we develop a general theory of topolog- 
ical defects and their associated gapless modes in Bloch 
and Bloch-BdG theories in all symmetry classes. As in 
Ref. [5^, we assume that far away from the defect the 
Hamiltonian varies slowly in real space, allowing us to 
consider adiabatic changes of the Hamiltonian as a func- 
tion of the real space position r. We thus seek to classify 
Hamiltonians 7^(14, r), where k is defined in a c? dimen- 
sional Brillouin zone (a torus T'*), and r is defined on 
a D dimensional surface 5"^ surrounding the defect. A 
similar approach can be used to classify cyclic tempo- 
ral variations in the Hamiltonian, which define adiabatic 
pumping cycles. Hereafter we will drop the BdG sub- 
script on the Hamiltonian with the understanding that 
the symmetry class dictates whether it is a Bloch or BdG 
Hamiltonian. 

In Fig. [T]we illustrate the types of topological defects 
that can occur in d = 1, 2 or 3. For D = we regard 5° 
as two points ({ — 1,+!}). Our topological classification 
then classifies the difference of 'H(k, -1-1) and ■H(k, — 1). A 
non trivial difference corresponds to an interface between 



two topologically distinct phases. For D = 1 the one 
parameter families of Hamiltonians describe line defects 
in (i = 3 and point defects in = 2. For d = 1 it 
could correspond to an adiabatic temporal cycle i?(k, t). 
Similarly for D = 2, the two parameter family describes 
a point defect for d = 3 or an adiabatic cycle for a point 
defects in d = 2. 

Classifying the D parameter families of d dimensional 
Bloch-BdG Hamiltonians subject to symmetries leads to 
a generalization of the periodic table discussed above. 
The original table corresponds to D = 0. For D > 
we find that for a given symmetry class the topological 
classification (Z, Z2 or 0) depends only on 



D. 



(1.1) 



Thus, all line defects with S = 2 have the same topolog- 
ical classification, irrespective of d, as do point defects 
with (5 = 1 and pumping cycles with 5 = 0. Though 
the classifications depend only on 5, the formulas for the 
topological invariants depend on both d and D. 

This topological classification of H(k,r) suggests a 
generalization of the bulk-boundary correspondence that 
relates the topological class of the Hamiltonian character- 
izing the defect to the structure of the protected modes 
associated with the defect. This has a structure remi- 
niscent of a mathematical index theorem^ that relates a 
topological index to an analytical index that counts the 
number of zero modes^i^i^i^"— . In this paper we will 
not attempt to prove the index theorem. Rather, we will 
observe that the topological classes for H(k,r) coincide 
with the expected classes of gapless defect modes. In 
this regards the dependence of the classification on S in 
p.ip is to be expected. For example, a point defect at 
the end of a one dimensional system (5=1—0) has the 
same classification as a point defects in two dimensions 
((5 = 2 — 1) and three dimensions ((5 = 3 — 2). 

We will begin in section [H] by describing the general- 
ized periodic table. We will start with a review of the 
Altland Zirnbauer symmetry classes^ and a summary of 
the properties of the table. In Appendix |^ we will jus- 
tify this generalization of the table by introducing a set 
of mathematical mappings that relate Hamiltonians in 
different dimensions and different symmetry classes. In 
addition to establishing that the classifications depend 
only on (5 = d — D, these mappings allow other features 
of the table, already present for 13 = to be easily un- 
derstood, such as the pattern in which the classifications 
vary as a function of symmetry class as well as the Bott 
periodicity of the classes as a function of d. 

In section IIIII and IIVI we will outline the physical con- 
sequences of this theory by discussing a number of ex- 
amples of line and point defects in different symmetry 
classes and dimensions. The simplest example is that of 
a line defect in a 3D system with no symmetries. In sec- 
tion [TTr^ we will show that the presence of a ID chiral 
Dirac fermion mode (analogous to an integer quantum 
Hall edge state) on the defect is associated with an inte- 
ger topological invariant that may be interpreted as the 
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winding number of the "0" term that characterizes the 
magnetoelectric polarizabihtj*i^. This description uni- 
fies a number of methods for "engineering" chiral Dirac 
fermions, which will be described in several illustrative 
examples. 

Related topological invariants and illustrative exam- 
ples will be presented in Sections IIII BIHiTeI for line de- 
fects in other symmetry classes that are associated with 
gapless ID helical Dirac fermions, ID chiral Majorana 
fermions and ID helical Majorana fermions. In section 
lIVI wc will consider point defects in ID models with chiral 
symmetry such as the Jackiw Rebbi model^ or the Su, 
Schrieffer, Heeger model^, and in superconductors with- 
out chiral symmetry that exhibit Majorana bound states 
or Majorana doublets. These will also be related to the 
early work of Jackiw and Rossi^ on Majorana modes at 
point defects in a model with chiral symmetry. 

Finally, in Section |V] we will regard r as including a 
temporal variable, and apply the considerations in this 
paper to classify cyclic pumping processes. The Thouless 
charge pump^i^ corresponds to a non trivial cycle in a 
system with no symmetries and S = {d = D = 1). A 
similar pumping scenario can be applied to superconduc- 
tors and defines a fermion parity pump. This, in turn, is 
related to the non-Abelian statistics of Ising anyons, and 
provides a framework for understanding braidless oper- 
ations on systems of three dimensional superconductors 
hosting Majorana fermion bound states. Details of sev- 
eral technical calculations can be found in the Appen- 
dices. 

An interesting recent preprint by Freedman et alj^, 
which appeared when this manuscript was in its final 
stages discusses some aspects of the classification of topo- 
logical defects in connection with a rigorous theory of 
non-Abelian statistics in higher dimensions. 



II. PERIODIC TABLE FOR DEFECT 
CLASSIFICATION 

Tabic H] shows the generalized periodic table for the 
classification of topological defects in insulators and su- 
perconductors. It describes the equivalence classes of 
Hamiltonians 7^(14, r), that can be continuously deformed 
into one another without closing the energy gap, subject 
to constraints of particle- hole and/or time reversal sym- 
metry. These arc mappings from a base space defined by 
(k, r) to a classifying space, which characterizes the set of 
gapped Hamiltonians. In order to explain the table, we 
need to describe (i) the symmetry classes, (ii) the base 
space, (iii) the classifying space and (iv) the notion of 
stable equivalence. The repeating patterns in the table 
will be discussed in Section III CI Much of this section is 
a review of material in Refs. I35ll36l . What is new is the 
extension to £> > 0. 
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TABLE I: Periodic table for the classification of topological 
defects in insulators and superconductors. The rows cor- 
respond to the different Altland Zirnbauer (AZ) symmetry 
classes, while the columns distinguish different dimensionali- 
ties, which depend only on & = d — D. 

A. Symmetry Classes 

The presence or absence of time reversal symmetry, 
particle- hole symmetry and/or chiral symmetry define 
the 10 Altland Zirnbauer symmetry classes^. Time re- 
versal symmetry implies that 

■H(k,r) = e-H(-k,r)e-\ (2.1) 

where the anti unitary time reversal operator may be 
written 8 = e^^^^^^K. S"^ is the spin and K is complex 
conjugation. For spin 1/2 fermions, = —1, which 
leads to Kramers theorem. In the absence of a spin or- 
bit interaction, the extra invariancc of the Hamiltonian 
under rotations in spin space allows an additional time 
reversal operator 6' = to be defined, which satisfies 

e'2 = +1. 

Particle-hole symmetry is expressed by 

H(k,r) = -SH(-k,r)S-i, (2.2) 

where S is the anti unitary particle-hole operator. Fun- 
damentally, — +1. However, as was the case for Q, 
the absence of spin orbit interactions introduces an addi- 
tional particle hole symmetry, which can satisfy = — 1 . 

Finally, chiral symmetry is expressed by a unitary op- 
erator H, satisfying 

H(k,r) = -HH(k,r)n-i. (2.3) 

A theory with both particle-hole and time reversal sym- 
metry automatically has a chiral symmetry H = e^-^OS. 
The phase x can be chosen so that H^ = 1. 

Specifying 6^ = 0, ±1, = 0, ±1 and H^ = 0, 1 (here 
denotes the absence of symmetry) defines the 10 Alt- 
land Zirnbauer symmetry classes. They can be divided 
into two groups: 8 real classes that have anti unitary sym- 
metries O and or S plus 2 complex classes that do not 
have anti unitary symmetries. Altland and Zirnbauer's 
notation for these classes, which is based on Cartan's 
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FIG. 2: The 8 Real symmetry classes that involve the anti 
unitary symmetries G (time reversal) and/or H (particle-hole) 
are specified by the values of = ±1 and — ±1. They 
can be visualized on an eight hour "clock" . 



classification of symmetric spaces, is shown in the left 
hand part of Table ID 

To appreciate the mathematical structure of the 8 real 
symmetry classes it is helpful to picture them on an 8 
hour "clock", as shown in Fig. [21 The x and y axes of 
the clock represent the values of and 9^. The "time" 
on the clock can be represented by an integer s defined 
modulo 8. Kitae\«2£ used a slightly different notation to 
label the symmetry classes. In his formulation, class D is 
described by a real Clifford algebra with no constraints, 
and in the other classes Clifford algebra elements are con- 
strained to anticommute with q positive generators. The 
two formulations are related by s = g -I- 2 mod 8. The 
complex symmetry classes can similarly be indexed by an 
integer s defined modulo 2. For all classes, the presence 
of chiral symmetry is associated with odd s. 



B. Base space, Classifying space and stable 
equivalence 

The Hamiltonian is defined on a base space composed 
of momentum k, defined in a d dimensional Brillouin zone 
T'^ and real space degrees of freedom r in a sphere (or 
S^~^ X for an adiabatic cycle). The total base space 
is therefore T'^ x (or T<^ x S^-^ x S'^). As in Ref. H, 
we will simplify the topological classification by treating 
the base space as a sphere 5"*+^. The "strong" topo- 
logical invariants that characterize the sphere will also 
characterize T"^ x . However, there may be additional 
topological structure in T"^ x that is absent in S"^^^. 
These correspond to "weak" topological invariants. For 
D = these arise in layered structures. A weak topologi- 
cal insulator, for example, can be understood as a layered 
two dimensional topological insulator. There are similar 
layered quantum Hall states. For _D ^ 0, then there will 
also be a "weak" invariant if the Hamiltonian H(k, Tq) for 
fixed r = rg is topologically non trivial. As is the case 
for the classification of bulk phases D = 0, we expect 



that the topologically protected gapless defect modes are 
associated with the strong topological invariants. 

The set of Hamiltonians that preserve the energy 
gap separating positive and negative energy states can 
be simplified without losing any topological informa- 
tion. Consider the retraction of the original Hamiltonian 
'H(k, r) to a simpler Hamiltonian whose eigenvalue spec- 
trum is "flattened" , so that the positive and negative en- 
ergy states all have the same energy zLEq. The flattened 
Hamiltonian is then specified the set of all n eigenvectors 
(defining a U (n) matrix) modulo unitary rotations within 
the k conduction bands or the n ~ k valence bands. The 
flattened Hamiltonian can thus be identifled with a point 
in the Grassmanian manifold. 



Gn,k U{n)/U{k) X U{n - fc). 



(2.4) 



It is useful to broaden the notion of topological equiv- 
alence to allow for the presence of extra trivial energy 
bands. Two families of Hamiltonians are stably equiv- 
alent if they can be deformed into one another after 
adding an arbitrary number of trivial bands. Thus, triv- 
ial insulators with different numbers of core energy levels 
are stably equivalent. Stable equivalence can be imple- 
mented by considering an expanded classifying space that 
includes an inflnite number of extra conduction and va- 
lence bands, Co = U/U x U = IJfc^o ^oo.fc- 

With this notion of stable equivalence, the equivalence 
classes of Hamiltonians H(k,r) can be formally added 
and subtracted. The addition of two classes, denoted 
[Hi] + [^2] is formed by simply combining two indepen- 
dent Hamiltonians into a single Hamiltonian given by 
the matrix direct sum, [Hi (B 'H2]. Additive inverses are 
constructed through reversing conduction and valence 
bands, ['Hi]-[H2] = [^10-^2]- is guaranteed 

to the be trivial class [0]. Because of this property, the 
stable equivalence classes form an Abelian group, which 
is the key element of K theorji^"— . 

Symmetries impose constraints on the classifying 
space. For the symmetry classes with chiral symmetry, 
(|2.3p restricts n = 2k and the classifying space to a sub- 
set Ci = U{oo) C U/U X U. The anti unitary symmetries 
(|2.1l2.2p impose further constraints. At the special points 
where k and — k coincide, the allowed Hamiltonians are 
described by the 8 classifying spaces TZq of Real K theory. 



C. Properties of the periodic table 

For a given symmetry class s, the topological classifi- 
cation of defects is given by the set of stable equivalence 
classes of maps from the base space (k, r) S 5^+^^ to the 
classifying space, subject to the symmetry constraints. 
These form the K group, which we denote as Kc{s; D, d) 
for the complex symmetry classes and Kr{s; D, d) for the 
real symmetry classes. These are listed in Table H] 

Table H] exhibits many remarkable patterns. Many can 
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be understood from the following basic periodicities, 

Kf{s]D,d+l) = Kf{s-l;D,d), (2.5) 
Kw{s;D + l,d) = Kp{s + l;D,d). (2.6) 

Here s is understood to be defined modulo 2 for F = 
C and modulo 8 for F = M. We will establish these 
identities mathematically in Appendix]^ The basic idea 
is to start with some Hamiltonian in some symmetry class 
s and dimensionalities D and d. It is then possible to 
explicitly construct two new Hamiltonians in one higher 
dimension which have either (i) d — > d + 1 or (ii) D — >■ 
D+l. These new Hamiltonians belongs to new symmetry 
classes that are shifted by one "hour" on the symmetry 
clock and characterized by (i) s — >■ s + 1 or (ii) s ^ s — 1. 
We then go on to show that this construction defines 
a 1-1 correspondence between the equivalence classes of 
Hamiltonians with the new and old symmetry classes and 
dimensions, thereby establishing (j2.5p and (j2.6p . 

The periodicities (j2.5|) and (|2.6p have a number of con- 
sequences. The most important for our present purposes 
is they can be combined to give 

Kw{s;D + l,d+l)=Kvis;D,d). (2.7) 

This (1,1) periodicity shows that the dependence on the 
dimensions d and D only occurs via i5 = d — D. Thus the 
dependence of the classifications on D can be deduced 
from the table for Z) = 0. This is one of our central 
results. 

In addition, the periodicities (|2.5|2.6|) explain other 
features of the table that are already present for D = 0. 
In particular, the fact that s is defined modulo 2 (8) for 
the complex (real) classes leads directly to the Bott pe- 
riodicity of the dependence of the classifications on d: 

Kcis;D,d + 2) = Kcis;D,d), (2.8) 
KK{s;D,d + 8) = KMis;D,d). (2.9) 

Moreover, (|2.5I2.6|) shows that Ka{s; D,d) depends only 
on d — D — s. This explains the diagonal pattern in 
Table HI in which the dependence of the classification on 
d is repeated in successive symmetry classes. Thus, the 
entire table could be deduced from a single row. 

Equations (|2.5p and (|2.6p do not explain the pattern 
of classifications within a single row. Since this is a 
well studied math problem there are many routes to the 
answe r^^i^^'^^ . One approach is to notice that for d = 0, 
Kr{s, D, 0) is simply the D'th homotopy group of the ap- 
propriate classifying space which incorporates the sym- 
metry constraints. For example, for class BDI (s = 1, 
'E? = +1, = +1) the classifying space is the orthogo- 
nal group 0(oo). Then, K]fi{l,D,0) = tt 13(0(00)), which 
are well known. This implies 

K^is; D, d) = ^,,+,5-d-i(0(oo)). (2.10) 

Additional insight can be obtained by examining the 
interconnections between different elements of the table. 



Symmetry Topological classes ID Gapless Fermion modes 



A Z Chiral Dirac 

D Z Chiral Majorana 

Dili Z2 Helical Majorana 

All Z2 Helical Dirac 

C 2Z Chiral Dirac 



TABLE II: Symmetry classes that support topologically non 
trivial line defects and their associated protected gapless 
modes. 



For example, the structure within a column can be an- 
alyzed by considering the effect of "forgetting" symme- 
tries. Hamiltonians belonging to the real chiral (non chi- 
ral) classes are automatically in complex class AIII (A). 
There are therefore K group homomorphisms that send 
any real entries in table|T]to complex ones directly above. 
In particular, as detailed in Appendix|B]this distinguishes 
the Z and 2Z entries, which indicate the possible values 
of Chern numbers (or U(n) winding numbers) for even 
(or odd) S. In addition, the dimensional reduction argu- 
ments given in Refs. IioIItiI lead to a dimensional hierar- 
chy, which helps to explain the pattern within a single 
row as a function of d. 



III. LINE DEFECTS 

Line defects can occur at the edge of a 2D system 
((5 = 2 - 0) or in a 3D system {6 = 3- 1). From TableU 
it can be seen that there are five symmetry classes which 
can host non trivial line defects. These are expected to 
be associated with gapless fermion modes bound to the 
defect. Table |lT] lists non trivial classes, along with the 
character of the associated gapless modes. In the follow- 
ing subsections we will discuss each of these cases, along 
with physical examples. 



A. Class A: Chiral Dirac Fermion 

1. Topological Invariant 

A line defect in a generic 3D Bloch band theory with 
no symmetries is associated with an integer topological 
invariant. This determines the number of chiral Dirac 
fermion modes associated with the defect. Since ■H(k, r) 
is defined on a compact 4 dimensional space, this invari- 
ant is naturally expressed as a second Chern number, 

n=^[ Tr[^A^], (3.1) 

where 

T = dA + AAA (3.2) 
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is the curvature form associated with the non-Abehan 
Berry's connection Aij = {ui\duj) characterizing the va- 
lence band eigenstates \uj{'k,s)) defined on the loop 
parameterized by s. 

It is instructive to rewrite this as an integral over s of 
a quantity associated with the local band structure. To 
this end, it is useful to write Tr[J^A J^] = dQa, where the 
Chern Simons 3 form is, 

Q3 ^ Ti-[A A dA+'^A /\ A /\ A]. (3.3) 

Now divide the integration volume into thin slices, x 
AS*^, where AS^ is the interval between s and s + As. 
In each slice. Stokes' theorem may be used to write the 
integral as a surface integral over the surfaces of the slice 
at s and s + As. In this manner, Eq. I3.1l mav be written 

Ztt Jgi as 

where 

^(■'*) = ^^3 23(k,s). (3.5) 

Eq. 13.51 is precisely the Qi, Hughes, Zhang formulaF^ for 
the "0" term that characterizes the magnetoelectric re- 
sponse of a band insulator. 9 = for an ordinary time 
reversal invariant insulator, and = tt in a strong topo- 
logical insulator. If parity and time reversal symmetry 
are broken then 6 can have any intermediate value. We 
thus conclude that the topological invariant associated 
with a line defect, which determines the number of chiral 
fermion branches is given by the winding number of 9. 

We now consider several examples of 3D line defects 
that are associated with chiral Dirac fermions. 



2. Dislocation in a 3D Integer quantum Hall state 

A three dimensional integer quantum Hall state can 
be thought of as a layered version of the two dimensional 
integer quantum Hall state. This can be understood most 
simply by considering the extreme limit where the layers 
are completely decoupled 2D systems. A line dislocation, 
as shown in Fig. [3] will then involve an edge of one of the 
planes and be associated with a chiral fermion edge state. 
Clearly, the chiral fermion mode will remain when the 
layers are coupled, provided the bulk gap remains finite. 
Here we wish to show how the topological invariant p.ip 
reflects this fact. 

On a loop surrounding the dislocation parameterized 
by s € [0,1] we may consider a family of Hamiltonians 
-ff (k, s) given by the Hamiltonian of the original bulk 
crystal displaced by a distance sB, where B is a lattice 
vector equal to the Burgers vector of the defect. The 
corresponding Bloch wavefunctions will thus be given by. 




FIG. 3: A line dislocation in a three dimensional quantum 
Hall state characterized by Burgers vector B. 

where Mmk(r) are Bloch functions for the original crystal. 
It then follows that the Berry's connection is 

A = A° + B-{k-a.P{k))ds, (3.7) 

where ^^„(k) = (w°„klVk|w°k> ' dk and 

a?5,„(k) = «kl(Vr + k)|<k). (3.8) 

With this definition, aP(k) is a periodic function: aP(k + 
G) = aP(k) for any reciprocal lattice vector GZS. 

If the crystal is in a three dimensional quantum Hall 
state, then the non zero first Chern number is an obstruc- 
tion to finding the globally continuous gauge necessary 
to evaluate p.5p . We therefore use p.ip . which can be 
evaluated by noting that 

Ti-[TAT] = Tr [B • {2T° A dk - d[T°, a^]) A ds] . (3.9) 

Upon integrating Tr[J^A J-"] the total derivative term van- 
ishes due to the periodicity of aF. Evaluating the integral 
is then straightforward. The integral over s trivially gives 
1. We are then left with 

n=-^B-G„ (3.10) 

ZTT 

where 

Gc = / dkATr[J-"]. (3.11) 

Gc is a reciprocal lattice vector that corresponds to the 
triad of Chern numbers that characterize a 3D system. 
For instance, in a cubic system Gc = {2t: / a)(nx,ny^nz), 
where, for example Uz = (27r)^^ j Tv[F^y\dkx A dky, for 
any value of kz- 

An equivalent formulation is to characterize the dis- 
placed crystal in terms of 9. Though p.5p can not be 
used, (|3.ip and p.4p can be used to implicitly define 9 
up to an arbitrary additive constant. 



(3.6) 



9[s) = sB • Gc. 



(3.12) 
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FIG. 4: Heterostructure geometries for chiral Dirac fermions. 
(a) and (b) show antiferromagnetic or ferromagnetic insula- 
tors on the surface of a topological insulator, with chiral Dirac 
fermions at a domain wall, (c) and (d) show a domain wall 
in an antiferromagnetic or ferromagnetic topological insula- 
tor. Chiral fermion modes are present when the domain wall 
intersects the surface. 



3. Topological insulator hetero structures 

Another method for engineering chiral Dirac fermions 
is use heterostructures that combine topological insula- 
tors and magnetic materials. The simplest version is a 
topological insulator coated with a magnetic film that 
opens a time reversal symmetry breaking energy gap at 
the surface. A domain wall is then associated with a 
chiral fermion mode. In this section we will show how 
this structure, along with some variants on the theme, 
fits into our general framework. We first describe the 
structures qualitatively, and then analyze a model that 
describes them. 

Fig. |4] shows four possible configurations. Figs. IHa,b) 
involve a topological insulator with magnetic materials 
on the surface. The magnetic material could be either 
ferromagnetic or antiferromagnetic. We distinguish these 
two cases based on whether inversion symmetry is broken 
or not. Ferromagnetism does not violate inversion sym- 
metry, while antiferromagnetism does (at least for inver- 
sion about the middle of a bond). This is relevant because 
, discussed above, is quantized unless both time reversal 
and inversion symmetries are violated. Of course, for a 
non centrosymmetric crystal inversion is already broken, 
so the distinction is unnecessary. 

Fig. lU^a) shows a topological insulator capped with 
antiferromagnetic insulators with 9 = ±€ separated by 
a domain wall. Around the junction where the three 
regions meet 9 cycles between tt, -\-e and — e. Of course 
this interface structure falls outside the adiabatic regime 
that p.5p is based on. However, it is natural to expect 



that the physics would not change if the interface was 
"smoothed out" with 9 taking the shortest smooth path 
connecting its values on either side of the interface. 

Fig. HJb) shows a similar device with ferromagnetic 
insulators, for which = or tt. In this case the adiabatic 
assumption again breaks down, however, as emphasized 
in Ref. [lol . the appropriate way to think about the surface 
is that 9 connects and tt along a path that is determined 
by the sign of the induced gap, which in turn is related to 
the magnetization. In this sense, 9 cycles by 27r around 
the junction. 

In Figs. |4l^c,d) we consider topological insulators which 
have a weak magnetic instability. If in addition to time 
reversal, inversion symmetry is broken, then ~ tt ± e. 
Recently Li, Wang, Qi and Zhang2^ have considered such 
materials in connection with a theory of a dynamical ax- 
ion and suggested that certain magnetically doped topo- 
logical insulators may exhibit this behavior. They re- 
ferred to such materials as topological magnetic insula- 
tors. We prefer to call them magnetic topological insula- 
tors because as magnetic insulators they are topologically 
trivial. Rather, they are topological insulators to which 
magnetism is added. Irrespective of the name, such ma- 
terials would be extremely interesting to study, and as we 
discuss below, may have important technological utility. 

Fig. Sfc) shows two antiferromagnetic topological in- 
sulators with = TT ± e separated by a domain wall, 
and Fig. Sfd) shows a similar device with ferromagnetic 
topological insulators. They form an interface with an 
insulator, which could be vacuum. Under the same con- 
tinuity assumptions as above the junction where the do- 
main wall meets the surface will be associated with a 
chiral fermion mode. Like the structure in Fig. |4|a), this 
may be interpreted as an edge state on a domain wall 
between the "half quantized" quantum Hall states of the 
topological insulator surfaces. However, an equally valid 
interpretation is that the domain wall itself forms a sin- 
gle two dimensional integer quantum Hall state with an 
edge state. Our framework for topologically classifying 
the line defects underlies the equivalence between these 
two points of view. 

Mong, Essen and Moore^^ have introduced a differ- 
ent kind of antferromagnetic topological insulator that 
relies on the symmetry of time reversal combined with 
a lattice translation. Due to the necessity of translation 
symmetry, however, such a phase is not robust to disor- 
der. They found that chiral Dirac modes occur at certain 
step edges in such crystals. These chiral modes can also 
be understood in terms of the invariant (|3.ip . Note that 
these chiral modes survives in the presence of disorder 
even though the bulk state does not. Thus, the chiral 
mode, protected by the strong invariant p.ip . is more 
robust than the bulk state that gave rise to it. 

If one imagines weakening the coupling between the 
two antiferromagnetic topological insulators (using our 
terminology, not that of Mong, et ali2^) and taking them 
apart, then at some point the chiral mode has to dis- 
appear. At that point, rather than taking the "shortest 
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path" between 7r±e, 9 takes a path that passes through 0. 
At the transition between the "short path" and the "long 
path" regimes, the gap on the domain wall must go to 
zero, allowing the chiral mode to escape. This will have 
the character of a plateau transition in the 2D integer 
quantum Hall effect. 

Structures involving magnetic topological insulators 
would be extremely interesting to study because with 
them it is possible to create chiral fermion states with a 
single material. Indeed, one can imagine scenarios where 
a magnetic memory, encoded in magnetic domains, could 
be read by measuring the electrical transport in the do- 
main wall chiral fermions. 

To model the chiral fermions in these structures we 
begin with the simple three dimensional model for trivial 
and topological insulators considered in Ref. 

■Ho = wMxa-k+ (TO + e|k|2)/i^. (3.13) 

Here S represents spin, and /i^ describes an orbital degree 
of freedom, m > describes the trivial insulator and m < 
describes the topological insulator. An interface where 
m changes sign is then associated with gapless surface 
states. 

Next consider time reversal symmetry breaking per- 
turbations, which could arise from exchange fields due to 
the presence of magnetic order. Two possibilities include 



^a/ — hafHy, 



n, 



(3.14) 
(3.15) 



Either haf or hf^z will introduce a gap in the surface 
states, but they have different physical content. Hp has 
an inversion symmetry given by 'Ho(k) = P'Ho(— k)P 
with P = fiz. Clearly, Ti / respects this inversion sym- 
metry. Haf does not respect P, but does respect PQ. 
We therefore associate Ji / with ferromagnetic order and 
Haf with antiferromagnetic order. 

Within the adiabatic approximation, the topological 
invariant p.ip can be evaluated in the presence of ei- 
ther p.l4p or p.lSp . The antiferromagnetic perturba- 
tion p.l4p is most straightforward to analyze because 
"Ho + Haf is a combination of 5 anticommuting Dirac 
matrices. On a circle surrounding the junction parame- 
terized by s it can be written in the general form 



i/(k,.s) =h(k,s)-7. 



(3.16) 



where 7 = {p-x<Jx, fJ-xCTy, (J-xCTz, fJ-z, fJ-y) and h(k, s) = 
(wk, m(s) -t- e|kp, haf{s)). For a model of this form, the 
second Chern number (|3.ip is given simply by the wind- 
ing number of the unit vector d(k, s) = h/|h| G S"* as 
a function of k and s. This is most straightforward to 
evaluate in the limit e — > 0, where d is confined to the 
"equator" (di, (i2, da, 0, 0) everywhere except near k ^ 
and |k| > 1/e. The winding number is determined by the 
behavior at k ~ 0, and may be expressed by p. 41) with 9 
given by 

i9 + ^f^af 



Kf 



(3.17) 



We therefore expect a topological line defect to oc- 
cur at an intersection between planes where m and haf 
change sign. The chiral fermion mode associated with 
this defect can seen explicitly if we solve a simple linear 
model, m = fzZ, haf = fyU- This model, which has the 
form of a harmonic oscillator, is solved in Appendix [Cj 
and explicitly gives the chiral Dirac fermion mode with 
dispersion 

E{kx)=VSgIi{fzfy)kx. (3.18) 



B. Class D: chiral Majorana fermions 

1. Topological invariant 

A line defect in a superconductor without time rever- 
sal symmetry is characterized by an integer topological 
invariant that determines the number of associated chi- 
ral Majorana fermion modes. Since the BdG Hamilto- 
nian characterizing a superconductor has the same struc- 
ture as the Bloch Hamiltonian, we can analyze the prob- 
lem by "forgetting" about the particle hole symmetry 
and treating the BdG Hamiltonian as if it was a Bloch 
Hamiltonian. The second Chern number, given by (|3.ip 
can be defined. It can be verified that any value of the 
Chern number is even under particle-hole symmetry, so 
that particle-hole symmetry does not rule out a non zero 
Chern number. We may follow the same steps as p.m3.5p 
to express the integer topological invariant as 



/ Tr[^A^], 



(3.19) 



where J- is the curvature form characterizing the BdG 
theory. As in p.4p . fi may be expressed as a winding 
number of 0, which is expressed as an integral over the 
Brillouin zone of the Chern Simons 3 form. The dif- 
ference between n and n is that fi characterizes a BdG 
Hamiltonian. If we considered the BdG Hamiltonian for 
a non superconducting insulator, then due to the dou- 
bling in the BdG equation, we would find 



2n. 



(3.20) 



In this case, the chiral Dirac fermion that occurs for a 27r 
[n = 1) winding of 9 corresponds to a 47r (n = 2) winding 
of 9. Superconductivity allows for the possibility of a 27r 
winding in 0: a chiral Dirac fermion can be split into a 
pair of chiral Majorana fermions. 



2. Dislocation in a layered topological superconductor 

The simplest example to consider is a dislocation in 
a three dimensional superconductor. The discussion 
closely parallels Section UlI A 21 and we find 



n = -!-B • Gc 
27r 



(3.21) 
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FIG. 5: Heterostructure geometries for Cliiral Majorana 
fermions. (a) shows a magnetic domain wall on the surface 
of a topological superconductor, while (b) shows an interface 
between a superconductor and a magnet on the surface of a 
topological insulator. 



where B is the Burgers veetor of the dislocation and 
characterizes the triad of first Chern numbers character- 
izing the 3D BdG Hamiltonian. A 3D system consist- 
ing of layers of a 2D topological superconductor will be 
characterized by a non zero Gc- Since, as a 3D super- 
conductor, the layered structure is in the topologically 
trivial class, such a state could be referred to as a weak 
topological superconductor. 

The simplest model system in this class is a stack of 2D 
Px + ipy superconductors. A dislocation would then have 
n = 1 and a single chiral Majorana fermion branch. A 
possible physical realization of the weak topological su- 
perconductor state is Sr2Ru04, which may exhibit triplet 
Px + ipy pairing. Since the spin up and spin down elec- 
trons make two copies of the spinless state, a dislocation 
will be associated with n = 2. Thus, we predict that 
there will be two chiral Majorana modes bound to the 
dislocation, which is the same as a single chiral Dirac 
fermion mode. 



3. Superconductor Heterostructures 

We now consider heterostructures with associated chi- 
ral Majorana modes. The simplest to consider is a BdG 
analog of the structures considered in Fig. [5] These 
would involve, for example an interface between a 3D 
time reversal invariant topological superconductor with 
a magnetic material with a magnetic domain wall. The 
analysis of such a structure is similar to that in Eq. 13.131 
if we replace the Pauli matrices describing the orbital de- 
gree of freedom fl with Pauli matrices describing Nambu 
space f . Protected chiral Majorana fermion modes of this 
sort on the surface of '^He-B with a magnetic domain wall 
have been recently discussed by Volovik^. 

In Ref. [13 a different method for engineering chiral 
Majorana fermions was introduced by combining an in- 
terface between superconducting and magnetic regions 
on the surface of a topological insulator. To describe this 
requires the 8 band model introduced in Ref. [ssl . 

H = Tzfix(J ■k + {m + e|kp)T2^2 -|- + hfXy. (3.22) 

(Here, for simplicity we consider only the antiferromag- 
netic term). The surface of the topological insulator oc- 
curs at a domain wall (say, in the x-y plane) where m(z) 



changes sign. The superconducting order parameter A 
and magnetic perturbation h both lead to an energy gap 
in the surface states. This Hamiltonian is straightforward 
to analyze because [H, Tj-/ij^] =0, which allows the 8x8 
problem to be divided into two 4x4 problems, which 
have superconducting/magnetic mass terms A ± /i. Near 
a defect where A = h the A + h gap never closes, while 
the A~h gap can be critical. A > h leads to a supercon- 
ducting state, while A < h leads to a quantum Hall like 
state. There is a transition between the two aX A = h. 

An explicit model for the line defect can be formulated 
with m{z) = fzZ, A — h = fyy and A+h = M. The topo- 
logical invariant p.l9[) can be evaluated using a method 
similar to p.l6p . and the chiral Majorana states can be 
explicitly solved along the lines of p.l8p . 



C. Class All: Helical Dirac fermions 

1. Topological Invariant 

Line defects for class AH are characterized by a Z2 
topological invariant. To develop a formula for this in- 
variant we follow the approach used in Ref. [T^to describe 
the invariant characterizing the quantum spin Hall insu- 
lator. 

As in the previous section, a line defect in three dimen- 
sions is associated with a four parameter space (k, r) € 
X S^. Due to time reversal symmetry, the second 
Chern number that characterized the line defects in p.ip 
must be zero. Thus there is no obstruction to defin- 
ing Bloch basis functions |u(k, r)) continuously over the 
entire base space. However, the time reversal relation 
between (— k, r) and (k, r) allows for an additional con- 
straint, so that the state is specified by the degrees of 
freedom in half the Brillouin zone. 

As in Ref. \7B it is useful to define a matrix 



i(k,r) = (u„(k, r)|e|u„(-k, r)). 



(3.23) 



Because |M„i(k, r)) and |u,i(— k, r)) are related by time 
reversal symmetry w(k, r) is a unitary matrix, that de- 
pends on the gauge choice for the basis functions. Locally 
it is possible to choose a basis in which 



w(k, r) ~ Wo, 



(3.24) 



where wq is independent of k and r, so that states at 
(±k, r) have a fixed relation. Since for k = w = ~uF , 
Wo must be antisymmetric. A natural choice is thus wo = 
i(j2 ® 1. 

The Z2 topological invariant is an obstruction to find- 
ing such a constrained basis globally. The constrained 
basis can be defined on two patches, but the basis func- 
tions on the two patches are necessarily related by a topo- 
logically non trivial transition function. In this sense, the 
Z2 invariant resembles the second Chern number in p.ip . 

In Appendix IeI we will generalize the argument de- 
veloped in Ref. |75| to show that the transition function 
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relating the two patches defines the Z2 topological invari- 
ant, which may be written^ 



All Tt[TAT]- f Q3 



mod 2, 



(3.25) 

where J- and Q3 arc expressed in terms of the Berry's 
connection A using p.2p and (|3.3p . The integral is over 
half of the base space (1/2)(T'^ x -S*^), defined such that 
(k, r) and (— k, r) are never both included. The second 
term is over the boundary of (l/2)(r^ x S^), which is 
closed under (k, r) (— k, r). Eq. 13.251 must be used 
with care because the Chern Simons form in the second 
term depends on the gauge. A different continuous gauge 
can give a different i/, but due to (|3.24[) . they must be 
related by an even integer. Thus, an odd number is dis- 
tinct. 

In addition to satisfying p.24p . it is essential to use 
a gauge in which at least Q3 is continuous on d^T^ x 
(though not necessarily on aU of ^T"^ x S^). This 
continuous gauge can always be found if the base space is 
a sphere 5^. However for T^xS^, the "weak" topological 
invariants can pose an obstruction to finding a continuous 
gauge. We will show how to work around this difficulty 
at the end of the following section. 



2. Dislocation in a weak topological insulator 

Ran, Zhang and Vishwanath recently studied the prob- 
lem of a line dislocation in a topological insulator—. 
They found that an insulator with non trivial weak topo- 
logical invariants can exhibit topologically protected he- 
lical modes at an appropriate line dislocation. In this 
section we will show that these protected modes are as- 
sociated with a non trivial Z2 invariant in p.25p . In 
addition to providing an explicit example for this in- 
variant, this formulation provides additional insight into 
why protected modes can exist in a weak topological in- 
sulator. As argued in Ref. lolfl^ the weak topological 
invariants lose their meaning in the presence of disor- 
der. The present considerations show that the helical 
modes associated with the dislocation are protected by 
the strong topological invariant associated with the line 
defect. Thus if we start with a perfect crystal and add 
disorder, then the helical modes remain, even though the 
crystal is no longer a weak topological insulator. The 
helical modes remain even if the disorder destroys the 
crystaline order, so that dislocations become ill defined, 
provided the mobility gap remains finite in the bulk crys- 
tal. In this case, the Hamiltonian has a non trivial wind- 
ing around the line defect, even though the defect has 
no obvious structural origin. Thus, the weak topological 
insulator provides a route to realizing the topologically 
protected line defect. But once present, the line defect is 
more robust than the weak topological insulator. 

To evaluate the Z2 invariant p.25p for a line disloca- 
tion we repeat the analysis in section IIII A 21 Because of 



the subtlety with the application of (|3.25|) we will first 
consider the simplest case of a dislocation in a weak topo- 
logical insulator. Afterwards we will discuss the case of 
a crystal with both weak and strong invariants. 

The Bloch functions on a circle surrounding a dis- 
location are described by (|3.6p . and the evaluation of 
Tr[J^ A J-] proceeds exactly as in p.7ll3.9|) . To evaluate 
the second term in p.25|) we need the Chern Simons 3 
form. One approach is to use p.3p and p.7p . However, 
this is not continuously defined on d{l/2){T^ x S^) be- 
cause A has a term B • kds that is discontinuous at the 
Brillouin zone boundary. An alternative is to write 



Q3 = Tr[B • {2A" A dk - [J"", a^]) A ds] 



(3.26) 



From p.9p this clearly satisfies Tr[J^ A J-] ~ dQ^, and it 
is defined continuously on 9(1/2) (T'^ x S^) as long as AP 
is continuously defined on 9(1/2)T^. For a weak topolog- 
ical insulator this is always possible, provided (l/2)r'^ is 
defined appropriately. Eq. 13.261 differs from Eq. I3.3l bv a 
total derivative. 

Combining ([53]) . ((X^ and ((X^ . the terms involv- 
ing aP cancel because aP is globally defined. (Note that 

is unchanged by a k dependent - but r independent ~ 
gauge transformation). This can not be said of the term 
involving A^ , however, because in a weak topological in- 
sulator is not globally defined on (1/2)T'^. Performing 
the trivial integral over s we then find 



where 



G, 



2lT 



TrfJ-"! A dk - 



-B • G„ mod 2, 



(3.27) 



r3 



Ty[A°] a dk. (3.28) 



The simplest case to consider is a weak topological 
insulator consisting of decoupled layers of 2D quantum 
spin Hall insulator stacked with a lattice constant a in the 
z direction. In this case = J-^{kx,ky) is independent 
of kz, so the kz integral can be performed trivially. This 
leads to Gi, ~ (27r/a)i/z, where 



2^ 



Tr[^0] 



Tr[^0] 



(3.29) 



is the 2D Z2 topological invariant characterizing the in- 
dividual layers. 

Eq. 13.281 also applies to a more general 3D weak topo- 
logical insulator. A weak topological insulator is charac- 
terized by a triad of Z2 invariants that define a 
mod 2 reciprocal lattice vector— li^, 



(3.30) 



where bj arc primitive reciprocal lattice vectors cor- 
responding to primitive lattice vectors a^ (such that 
a^ • hj = 2Tr6ij). The indices Vi can be determined by 
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evaluating the 2D invariant p.29p on the time reversal 
invariant plane k • = tt. 

To show that G„ in (|3.28|) and (|3.30p are equivalent, 
consider G^, ■ ai in p.28p . If we write k Xihi + X2h2 + 
xabs, then the integrals over X2 and X3 have the form of 
(j3.29|) . Since this is quantized, it must be independent of 
Xi, and will be given by its value at Xi = 1/2. This then 
gives • ai = 2Tn^i. A similar analysis of the other 
components establishes the equivalence. A non trivial 
value of Eq. 13.271 is the same as the criterion for the 
existence of protected helical modes on a dislocation Ran, 
Zhang, Vishwanath^ derived using a different method. 

Evaluating p.28p in a crystal that is both a strong topo- 
logical insulator and a weak topological insulator (such as 
Bii-ajSbj;) is problematic because the 2D invariants eval- 
uated on the planes a;i = and xi = 1/2 are necessarily 
different in a strong topological insulator. This arises be- 
cause a non trivial strong topological invariant j/q is an 
obstruction to continuously defining on 9(1/2)T'^, so 
p.29|) can not be evaluated continuously between xi = 
and xi = 1/2. From the point of view of the topological 
classification of the defect on x S^, vq is like a weak 
topological invariant because it a property of and is 
independent of the real space parameter s in 5^. Thus 
this complication is a manifestation of the fact that topo- 
logical classification of Hamiltonians on T"^ x has more 
structure than those on S*^. The problem is not with the 
existence of the invariant i> on X S*!, but rather with 
applying the formulas p.25l3.28p . The problem can be 
circumvented with the following trick. 

Consider an auxiliary Hamiltonian H(k, r) = H(k, r)© 
UsTiO^), where JisTi is a simple model Hamiltonian 
for a strong topological insulator like Eq. I3.13[ which 
can be chosen such that it is a constant independent of 
k everywhere except in a small neighborhood close to 
k = where a band inversion occurs. Adding such a 
Hamiltonian that is independent of r will have no effect 
on the topologically protected modes associated with a 
line defect, so we expect the invariant v to be the same 
for both H(k,r) and H(k,r). If 7^(k, r) has a non triv- 
ial strong topological invariant i/q = 1 then ■H(k, r) will 
have vq = 0, so that Eq. 13.281 can be applied. Gi, will 
then be given by the 2D invariant p.29p evaluated for 
H, which will be independent of xi. Since JisTiO^) is k 
independent everywhere except a neighborhood of k = 0, 
this will agree with the 2D invariant evaluated for H at 
xi = 1/2, but not xi = 0. It then follows that even in a 
strong topological insulator the invariant characterizing 
a line dislocation is given by (|3.27p . where is given 
by p.30p in terms of the weak topological invariants. 



3. Hetero structure geometries 

In principle it may be possible to realize ID helical 
fermions in a 3D system that does not rely on a weak 
topological insulating state. It is possible to write down 
a 3D model, analogous to p.l3p that has bound helical 




FIG. 6: Helical Majorana fermions at a linear Josephson junc- 
tion with phase difference tt on the surface of a topological 
insulator. 



modes. However, it is not clear how to physically im- 
plement this model. This model will appear in a more 
physical context as a BdG theory in the following section. 



D. Class Dili: helical Majorana fermions 

Line defects for class Dili are characterized by a Z2 
topological invariant that signals the presence or absence 
of ID helical Majorana fermion modes. As in Section 
IIIIB[ the BdG Hamiltonian has the same structure as a 
Bloch Hamiltonian, and the Z2 invariant can be deduced 
by "forgetting" the particle hole symmetry, and treating 
the problem as if it was a Bloch Hamiltonian in class AH. 

There are several ways to realize helical Majorana 
fermions. The simplest is to consider the edge of a 2D 
time reversal invariant superconductor or superfluid, or 
equivalently a dislocation in a layered version of that 2D 
state. A second is to consider a topological line defect 
in a 3D class Dili topological superconductor or super- 
fluid. Such line defects are well known in of '^He B^i^ 
and have recently been revisited in Ref. [37ll43l . 

Here we will consider a different realization that uses 
topological insulators and superconductors. Consider a 
linear junction between two superconductors on the sur- 
face of a topological insulator as shown in Fig. [HI In Ref. 
[50I it was shown that when the phase difference between 
the superconductors is tt there are gapless helical Majo- 
rana modes that propagate along the junction. This can 
be described by an 8 band minimal model that describes 
a topological insulator surface with a superconducting 
proximity effect. 



U = vT2fi.j;a ■ k + {m + e|k|^)r^^^ 4- Ait^ 



(3.31) 



Here m is the mass describing the band inversion of a 
topological insulator, as in p. 221) . and Ai is the real part 
of the superconducting gap parameter. This model has 
time reversal symmetry with Q = iUyK and particle- 
hole symmetry with S = UyTyK. The imaginary part 
of the superconducting gap, A.2Ty violates time reversal 
symmetry. A line junction along the x direction with 
phase difference tt at the surface of a topological insulator 
corresponds to the intersection of planes where m(z) and 
Ai(y) change sign. 

The Z2 invariant characterizing such a line defect is 
straightforward to evaluate because [■H,/ij,Ta;] = 0. This 
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extra symmetry allows a "spin Chern number" to be de- 
fined, ria- = (167r^)~^ / Tr[/ZyTa; J" A J"]. Since the sys- 
tem decouples into two time reversed versions of p.l3p . 
rig- = 1. By repeating the formulation in Appendix lEl of 
the Z2 invariant i/, it is straightforward to show that this 
means 1^=1. 

The helical modes can be explicitly seen by solving the 
linear theory, m = f^z, Ai = /j,y, which leads to the 
harmonic oscillator model studied in Appendix [C] In 
the space of the two zero modes the Hamiltonian has the 
form 

n^vk^a.,, (3.32) 
and describes ID helical Majorana fermions. 

E. Class C: Chiral Dirac fermions 

We finally briefly consider line defects in class C. Class 
C can be realized when time reversal symmetry is bro- 
ken in a superconductor without spin orbit interactions 
that has even parity singlet pairing. Line defects are 
characterized by an integer topological invariant that de- 
termines the number of chiral Majorana fermion modes 
associated with the line. As in class D, this may be evalu- 
ated by "forgetting" the particle-hole symmetry and eval- 
uating the corresponding Chern number that would char- 
acterize class A. The 2Z in Table[n]for this case, however, 
means that the Chern integer computed in this manner is 
necessarily even. This means that there will necessarily 
be an even number 2n of chiral Majorana fermion modes, 
which may cquivalcntly viewed as n chiral Dirac fermion 
modes. 

An example of such a system would be a 2D 
dx'i-y'i + idxy superconductor—, which exhibits chiral 
Dirac fermion edge states, or equivalently a dislocation 
in a 3D layered version of that state. 



IV. POINT DEFECTS 

Point defects can occur at the end of a ID system {S = 
1 — 0) or at topological defects in 2D (5 = 2—1) or 3D 
{6 — 3 — 2) systems. From the 6=1 column Table [H 
it can be seen that there arc five symmetry classes that 
can have topologically non trivial point defects. These 
are expected to be associated with protected zero energy 
bound states. Table Hill lists the non trivial classes, along 
with the character of the associated zero modes. In this 
section we will discuss each of these cases. 



A. Classes AIII, BDI and CII: chiral zero modes 

1. Topological invariant and zero modes 

Point defects in classes AIII, BDI and CII are charac- 
terized by integer topological invariants. The formula for 



Symmetry Topological classes E = Q Bound States 



AIII 
BDI 

D 
Dili 

CII 



Z 
Z 
Z2 
Z2 

2Z 



Chiral Dirac 
Chiral Majorana 
Majorana 

Majorana Kramers Doublet 
{= Dirac) 

Chiral Majorana Kramers 
Doublet (=Chiral Dirac) 



TABLE III: Symmetry classes supporting non trivial point 
topological defects and their associated E = Q modes. 



this integer invariant can be formulated by exploiting the 
chiral symmetry in each class. In a basis where the chiral 



symmetry operator is 11 
written. 



Tz, the Hamiltonian may be 



H(k,r) 





g(k,r)t 



9(k,r) 




(4.1) 



When the Hamiltonian has a flattened eigenvalue spec- 
trum Ti,^ = 1, g(k, r) is a unitary matrix. For a point 
defect in d dimensions, the Hamiltonian as a function of 
d momentum variables and D = d — 1 position variables 
is characterized by the winding number associated with 
the homotopy TT2d-i[U{n — > 00)] = Z, which is given by 



id -1)1 



{2d - l)l{2TTi)d Jrpa^ga- 



Tv[{qdq^) 



U2d-l-i 



(4.2) 



For a Hamiltonian that is built from anticommuting 
Dirac matrices, ■H(k, r) = d(k, r) • 7, this invariant 
is given simply by the winding degree of the mapping 
d(k, r) from T'^ x S'^~^ to S^'^~^, which is expressed as 
an integral of the Jacobian, 



2Tr<i 



(4.3) 



In class AIII there are no constraints on q{'k,r) other 
than unitarity, so all possible values of n are possible. 
There are additional constraints for the chiral classes 
with antiunitary symmetries. As shown in Appendix |B1 
this is simplest to see by analyzing the constraints on 
the winding degree discussed above, n must be zero in 
classes CI and DHL There is no constraint on n in class 
BDI, while n must be even in class CII. 

The topological invariant is related to an index that 
characterizes the chirality of the zero modes, 



(4.4) 



where N± are the number of zero modes that are eigen- 
states of H with eigenvalue ±1. To see that these zero 
modes are indeed protected consider A^-|_ = 71 > and 
A„ = 0. Any term in the Hamiltonian that could shift 
any of the N+ degenerate states would have to have a 
nonzero matrix element connecting states with the same 
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chirality. Such terms are forbidden, though, by the chi- 
ral symmetry {H, 11} = 0. In the superconducting classes 
BDI and CII the zero energy states are Majorana bound 
states. In class CII, however, since time reversal sym- 
metry requires that n must be even, the paired Majo- 
rana states can be regarded as zero energy Dirac fermion 
states. 

In the special case where H(k, r) has the form of a mas- 
sive Dirac Hamiltonian, by introducing a suitable regu- 
larization for |k| oo the topological invariant (|4.2l4.3p 
can be expressed in a simpler manner as a topological 
invariant characterizing the mass term. In the following 
subsections we consider this in the three specific cases 
1,2,3. 



fiTz- Combined with particle-hole symmetry ^ = UyTyK 
{"E? = +1), Q defines the BDI class with chiral symmetry 
n = r^. 

Evaluating the topological invariant (|4.2p again re- 
quires a |k| — >■ (X) regularization. One possibility is to 
add e|kpTa;, so that |k| — >■ oo can be replaced by a single 
point. In this case the invariant can be determined by 
computing the winding degree of d(k, r) on S^. In the 
limit e — ?> the k integral can be performed, so that (|4.2p 
can be expressed as the winding number of the phase of 



2. Solitons in d=l 

The simplest topological zero mode occurs in the 
Jackiw Rebbi model^, which is closely related to the 
Su Schrieffer Heeger model^. Consider 

x) = vkux + muy. (4-5) 

Domain walls where m{x) changes sign as a function of 
X are associated with the well known zero energy soliton 
states. 

To analyze the topological class requires a regulariza- 
tion for |fc| oo. This can cither be done with a lattice, 
as in the Su, Schrieffer, Heeger model or by adding a term 
efc^tTy, as in (|3.13p so that |fc| — >■ oo can be replaced by 
a single point. In either case, the invariant ()4.2|) changes 
by 1 when m changes sign. 

3. Jackiw Rossi Model m d=2 

Jackiw and Rossi introduced a two dimensional model 
that has protected zero modes^. The Hamiltonian can 
be written 

■H(k,r) =w7-k4-f-0(r), (4.6) 

where k = [kx, ky) and (71, 72) and (Fi, r2) are anticom- 
muting Dirac matrices. They showed that the core of a 
vortex where (f> = (f>i + i4>2 winds by 27rn is associated 
with n zero modes that are protected by the chiral sym- 
metry. Viewed as a BdG Hamiltonian, these zero modes 
are Majorana bound states. 

This can be interpreted as a Hamiltonian describing su- 
perconductivity in Dirac fermions. In this interpretation 
the Dirac matrices are expressed as (71,72) = Tz{ax,<yy) 
and (ri,r2) ~ {Tx,Ty), where (? is a Pauli matrix de- 
scribing spin and r describes particle-hole space. The 
superconducting pairing term is A = + i(j)2. In this 
interpretation a vortex violates the physical time reversal 
symmetry Q = iayK. However, even in the presence of a 
vortex this model has a fictitious "time reversal symme- 
try" Q = axTxK which satisfies 8^ = -1-1. This symme- 
try would be violated by a finite chemical potential term 



4- Hedgehogs in d = 3 

In Ref. [5^ we introduced a three dimensional model 
for Majorana bound states that can be interpreted as 
a theory of a vortex at the interface between a super- 
conductor and a topological insulator. In the special 
case that the chemical potential is equal to zero, model 
has the same form as (|4.6p . except that now all of the 
vectors are three dimensional. In the topological in- 
sulator model we have 7 = (71,72,73) = M^t^it and 
r = (ri,r2,r3) = {lJ,zTz,Tx,Ty). r and a arc defined 
as before, while jl describes a orbital degree of freedom. 
The chiral symmetry, II = fiyTz is violated if a chemical 
potential term fir^ is included. 

Following the same steps that led to (|4.7p the invariant 
(|4.2p is given by the winding number of = 0/101 on 5^, 

n= — [ • (d0 X d(f)). (4.8) 

B. Class D: Majorana bound states 

1. Topological invariant 

Point defects in class D are characterized by a Z2 topo- 
logical invariant that determines the presence or absence 
of a Majorana bound state associated with the defect. 
These include the well known end states in a ID p wave 
superconductor and vortex states in a 2d p^. -f ipy super- 
conductor. In Ref. [ill we considered such zero modes 
in a three dimensional BdG theory describing Majorana 
zero modes in topological insulator structures. Here we 
develop a unified description of all of these cases. 

For a point defect in d dimensions, the Hamiltonian 
depends on d momentum variables and D = d— 1 position 
variables. In appendix [Dl we show that the Z2 invariant 
is given by, 

^=lif^) / Q2d-imod2, (4.9) 
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where Q2d-i is the Chern Simons form. The specific 
cases of interest are, 

Si = Tr[^], (4.10) 
Qs - TiiAdA+'^A% (4.11) 

Qs - TiiAidAf + lA^dA+lA% (4.12) 
2 5 

It is instructive to see that (|4.9p reduces to (|4.2p in the 
case in which a system also has particle-hole symmetry. 
In this detailed in Appendix |D] it is possible to 

choose a gauge in which A = q^dq/2, so that Q2d-i oc 
(gdgt)2rf-i. 



^. End States m a ID superconductor 

The simplest example of a point defect in a supercon- 
ductor occurs in Kitaev's model^ of a one dimensional p 
wave superconductor. This is described by a simple ID 
tight binding model for spinless electrons, which includes 
a nearest neighbor hopping term tclci^i+h.c. and a near- 
est neighbor p wave pairing term AciCi+i -I- h.c. The 
Bogoliubov de Gennes Hamiltonian can then be written 
as 

n{k) = (tcosfc - /i)r^ + Asinfcr^. (4.13) 

This model exhibits a weak pairing phase for < t and 
a strong pairing phase for \fi\ > t. The weak pairing 
phase will have zero energy Majorana states at its ends. 

The topological invariant (|4.9p can be easily evalu- 
ated. We find A = d0/2, where 9 is the polar angle 
of d(fc) = {t cos k~ fi,A sin k). It follows that for < t, 
the topological invariant is v = 1 mod 2. 



3. Vortex in a 2D topological superconductor 

In two dimensions, a Majorana bound state occurs at a 
vortex in a topological superconductor. This can be eas- 
ily seen by considering the edge states of the topological 
superconductor in the presence of a hole^. Particle-hole 
symmetry requires that the quantized edge states come in 
pairs. When the flux is an odd multiple of h/2e, the edge 
states are quantized such that a zero mode is present. 
In this section we will evaluate the topological invariant 
(|4.9p associated with a loop surrounding the vortej*^^. 

We begin with the class D BdG Hamiltonian 
'Hp(kx,ky) characterizing the topological superconduc- 
tor when the superconducting phase is zero. We include 
the subscript p to denote the first Chern number that 
classifies the topological superconductor. We can then 
introduce a nonzero superconducting phase by a gauge 
transformation, 

Hp(k, ^) = e-*'^"^/2^"(k)e*'^"-/2, (4.14) 



where operates in the Nambu particle-hole space. We 
now wish to evaluate (j4.9p for this Hamiltonian when 
phase </3(s) winds around a vortex. There is, however, 
a difficulty because the Chern Simons formula requires 
a gauge that is continuous throughout the entire base 
space x S^. The nonzero Chern number p charac- 
terizing 'Hp(k) is an obstruction to constructing such a 
gauge. A similar problem arose in Section HHCH when 
we discussed a line dislocation in a weak topological su- 
perconductor. We can adapt the trick we used there to 
get around the present problem. We thus double the 
Hilbert space to include two copies of our Hamiltonian- 
one with Chern number p and one with Chern number 

^"W = (T\yk))- (4.15) 

We then put the vortex in only the -^p component, 

n{\i, if) = e-*'^«?i"(k)e''^'?, (4.16) 

where 

We added an extra phase factor by replacing by 1 -I- 
in order to make e"^' periodic under -|- 27r. 

Since the Chern number characterizing ■H°(k) is zero, 
there exists a continuous gauge, 

|ji.(k,^)) =e^^^|M°(k)), (4.18) 

which allows us to evaluate the Chern Simons integral. 
The Berry's connection Aij = {ui\duj) is given by 

A = A^ ^iQdip, (4.19) 

where ^°(k) is the connection describing 'H*'(k) and 
(5y(k) = (u°(k)|g|M°(k)). Inserting this into (|IT1|) and 
rearranging terms we find 

Qa = Tr[2gj-° - d{QA^)\ A dy^, (4.20) 

where = dA^^A^hA^. Since the second term is a to- 
tal derivative it can be discarded. For the first term there 
are two contributions from the 1 and the in (j4.17p . 
Upon integrating over k, the term can be shown to 
vanish as a consequence of particle-hole symmetry. The 
1 term simply projects out the Berry curvature of the 
original Hamiltonian Hp(k), so that 

Qs = Tr[J-°] A dip. (4.21) 

It follows from (|4.9p that the Z2 invariant characterizing 
the vortex is 

V = pm mod 2, (4.22) 

where p is the Chern number characterizing the topolog- 
ical superconductor and m is the phase winding number 
associated with the vortex. 
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FIG. 7: Visualization of the Z2 Pontrjagin invariant charac- 
terizing maps from (k, 0) G x to when the winding 
degree for each </> is 1. The inner sphere corresponds to k at 
<j> — 0, while the outer sphere is = 2n. The lines depict 
inverse images of two specific points on S^, which are lines 
connecting the inner and outer spheres. In (a) they have one 
twist, which can not be eliminated. The double twist in (b) 
can be unwound by smoothly dragging the paths around the 
sphere to arrive at (c), which has no twist. 

It is also instructive to consider this invariant in the 
context of the simple two band model introduced by Read 
and Greenes. This can be written as a simple tight bind- 
ing model 

T-L^lkxTky) = (t(cos fcj;+cos fcj,)— /i)T2+A(sin fc3;Ta;-|-sin fcj,Tj,). 

(4.23) 

where the superconducting order parameter A is real. 
As in (|4.13p , this model exhibits weak and strong pairing 
phases for < t and > t. These arc distinguished 
by the Chern invariant, which in turn is related to the 
winding number on of the unit vector d(k), where 

(i(k) are the coefficients of f in (|4.23p . A nonzero super- 
conducting phase is again introduced by rotating about 
Tz as in (|4.14l) . Here we wish to show that in this two 
band model the Z2 invariant v can be understood from 
a geometrical point of view. 

The Z2 invariant characterizing a vortex can be under- 
stood in terms of the topology of the maps d(kx, ky,(j)). 
from X to S^. These maps were first classified 
by Poutrjagin2^, and have also appeared in other phys- 
ical contexts^i^i^. Without losing generality, wc can 
reduce the torus to a sphere 5^, so the mappings are 
S"^ X ^ S'^. When for fixed cj) d{kx,ky,(f)) has an 
winding number of ±p, the topological classification 
is Z2p. In the case of interest, p = 1, so there arc two 
classes. 

This Z2 Pontrjagin invariant can be understood picto- 
rially by considering inverse image paths in (k, (f) space, 
which map to two specific points on S^. These corre- 
spond to ID curves in x S*^. Fig. [7] shows three 
examples of such curves. The inner sphere corresponds 
to = 0, while the outer sphere corresponds to (p = 2tt. 
Since p ~ 1, for every point on S" the inverse image path 
is a single curve connecting the inner and outer spheres. 
The key point is to examine the linking properties of 
these curves. The Z2 invariant describes the number of 
twists in a pair of inverse image paths, which is 1 in (a), 
2 in (b) and in (c). The configuration in (b) can be 
continuously deformed into that in (c) by dragging the 
paths around the inner sphere. This can be verified by a 



simple demonstration using your belt. The twist in (a), 
however, can not be undone. The number of twists thus 
defines the Z2 Pontrjagin invariant. 

4- Superconductor Hetero structures 

Finally, in three dimensions, a non trivial point de- 
fect can occur at a superconductor heterostructure. An 
example is a vortex in the superconducting state at the 
interface between a superconductor and a topological in- 
sulator. As shown in Ref. [5^ this can be described by 
the simple Hamiltonian, 

v. = vT^^^.d ■ k - + (m + e|kp)T2^2 + Ait^. + AiTy. 

(4.24) 

Here m is a mass which distinguishes a topological in- 
sulator from a trivial insulator, and A = Ai -t- zA2 is a 
superconducting order parameter. For p, — t)^ this Hamil- 
tonian has the form of the three dimensional version of 
(|4.6p discussed in IV. A. 4, where the mass term is char- 
acterized by the vector <\) = (m, Ai, A2). A vortex in A 
at the interface where m changes sign then corresponds 
to a hedgehog singularity in 0. From (|4.3|) . it can be 
seen that the class BDI Z invariant is ?i = 1. This then 
establishes that the class D Z invariant is = 1. The Z2 
survives when a nonzero chemical potential reduces the 
symmetry from class BDI to class D. 

C. Class Dili: Majorana doublets 

Point defects in class Dili are characterized by a Z2 
topological invariant. These are associated with zero 
modes, but unlike class D, the zero modes are required 
by Kramers theorem to be doubly degenerate. The zero 
modes thus form a Majorana doublet, which is equivalent 
to a single Dirac fermion. 

In Table m Class DHL 5 = 1 is an entry that is similar 
to Class All, (5 = 2. The Z2 for DIII invariant bears a 
resemblance to the invariant for All, which is a general- 
ization of the Z2 invariant characterizing the 2D quantum 
spin Hall insulator. In Appendix IE 21 we will establish a 
formula that employs the same gauge constraint, 

w(k,r)=wo, (4.25) 

where wq is a constant independent of k and r. i(j(k, r) 
relates the time reversed states at k and — k and is given 
by (|3.23p . Provided we choose a gauge that satisfies this 
constraint, the Z2 invariant is given by, 

^ = 7i (^) I 22d-i mod 2. (4.26) 

This formula is almost identical to the formula for a point 
defect in class D, but they differ by an important factor 
of two. Due to the combination of time reversal and 



particle-hole symmetry the Chern Simons integral (|4.26p 
is guaranteed to be an integer, but the integer is not 
gauge invariant. When the time reversal constraint is 
satisfied, the parity i> is gauge invariant. It then follows 
that the class D invariant in (|4.9|) . = mod 2. 

In the special case d = 1 there is a formula that does 
not rely on the gauge constraint, though it still requires a 
globally defined gauge. It is related to the similar "fixed 
point" formula for the invariant for the 2D quantum spin 
Hall insulator—, and has recently been employed by Qi, 
Hughes and Zhang^ to classify one dimensional time re- 
versal invariant superconductors. In class Dili, it is pos- 
sible to choose a basis in which the time reversal and par- 
ticle hole operators are given by 9 = TyK and S = TxK, 
so that the chiral operator is H = Tz. In this basis, the 
Hamiltonian has the form (|4.ip . where g(k, r) — )• q{k) sat- 
isfies q{—k) = —q{k)'^. Thus, Pf[(7(fc)] is defined for the 
time reversal invariant points k = and fc = tt. q(k) is 
related to 'w{k) because in a particular gauge it is possible 
to choose w{k) = q{k)/^y\Det [q{k)]\. The Z2 invariant is 
then given by, 

,^ Pf[g(7r)] VDct[g(0)] 
^ ^ Pf[g(0)] VDet[g(^)]' ^ ■ ' 

where the branch •y/Det[g(fc)] is chosen continuously be- 
tween fc = and k = IT. The equivalence of (|4.26p and 
(|4.27p for (i = 1 is demonstrated in Appendix IE 2 al Un- 
like (|4.26p . however, the fixed point formula (|4.27|) docs 
not have a natural generalization for d > 1. 

Majorana doublets can occur at topological defects in 
time reversal invariant topological superconductors, or in 
Helium 3B. Here we consider a different configuration at 
a Josephson junction at the edge of a quantum spin Hall 
insulator (Fig. |8]). When the phase difference across the 
Josephson junction is tt, it was shown in Rcf. Isslfssl that 
there is a level crossing in the Andreev bound states at 
the junction. This corresponds precisely to a Majorana 
doublet. 

This can be described by a the simple continuum ID 
theory introduced in Ref. |8^. 

n = vkT.a, + AiT^ (4.28) 

Here az describes the spin of the quantum spin Hall edge 
state, and Ai is the real superconducting order parame- 
ter. This model has particle-hole symmetry S = ciyTyK 
and time reversal symmetry Q = iayK and is in class 
DHL A TT junction corresponds to a domain wall where 
Al changes sign. Following appendix [Cl it is straight- 
forward to sec that this will involve a degenerate pair 
of zero modes indexed by the spin cr^ and chirality Ty 
constrained by Tya, = —1. 

The Hamiltonian (|4.28p should be viewed as an a low 
energy theory describing the edge of a 2D quantum spin 
Hall insulator. Nonetheless, we may describe a domain 
wall where Ai changes sign using an effective one dimen- 
sional theory by introducing a regularization replacing 
Al by Al -f- efc^. This regularization will not effect the 
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FIG. 8: A Josephson junction in proximity with the helical 
edge states of a quantum spin Hall insulator. When the phase 
difference is tt, there is a zero energy Majorana doublet at the 
junction. 

topological structure of a domain wall where Ai changes 
sign. A topologically equivalent lattice version of the 
theory then has the form. 

U = t sin kr^cT:, + (Ai + u{l - cosfc))r:i,. (4.29) 

where we assume |Ai| < 2u. 

The topological invariant can be evaluated using ei- 
ther (j4:26l) or (1127]). To use (gSSl), note that (ITM]) 
has exactly the same form as two copies (distinguished 
by cr^ = ±1) of ((4T31) . The evaluation of (|4:26)) then 
proceeds along the same lines. It is straightforward to 
check that in a basis where the time reversal constraint 
(|4.25p is satisfied (this fixes the relative phases of the 
az — ±1 states), A ~ d9, where 9 is the polar angle of 
d(fc) = (t sin fc, Al -|-m(1 — cos/c)). It follows that a defect 
where Ai changes sign has i> = 1. 

To use (j4.27|) . we transform to a basis in which 
= 'TyK, S = TxK and H = t^. This is 
accomplished by the unitary transformation U = 
exp[i(7r/4)f7j,T2] exp[i(7r/4)Tj;]. Then, T-l has the form 
of Eq. [5TT31 with q{k) = ~i{t shikar + (Ai -t- u(l - 
cosfc))crj,. It follows that dei[q{k)] is real and positive 
for aU k. Moreover, Pf [(?(0)]/A/det[g(0)] = sgn[Ai] while 
Pf [(7(7r)]/\/det[g(7r)] = 1. Again, a defect where Ai 
changes sign has D ~ 1. 



V. ADIABATIC PUMPS 

In this section we will consider time dependent Hamil- 
tonians H(k, r, t), where in additional to having adiabatic 
spatial variation r there is a cyclic adiabatic temporal 
variation parameterized by t. We will focus on point like 
spatial defects, in which the dimensions of k and r are 
related by d — Z? = 1 . 

Adiabatic cycles in which 'H(k, r,t ^ T) = 7i(k,r,t ^ 
0) can be classified topologically by considering t to be an 
additional "spacelike" variable, defining D = D + 1. Such 
cycles will be classified by the i5 = column of Table ID 
Topologically non trivial cycles correspond to adiabatic 
pumps. Table HVl shows the symmetry classes which host 
non trivial pumping cycles, along with the character of 
the adiabatic pump. There are two general cases. Classes 
A. AI and AH define a charge pump, where after one cycle 
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Symmetry Topological classes Adiabatic Pump 



A Z Charge 

AI Z Charge 

BDI Z2 Fermion Parity 

D Z2 Fermion Parity 

All 2Z Cliarge Kramers Doublet 



to the edge states surrounding the hole. This pumping 
process can be characterized by the second Chern number 
characterizing the 2D Hamiltonian fcy, 0, t), where 

9 parameterizes a circle surrounding the hole. A similar 
pump in 3D can be considered, and is characterized by 
the third Chcrn number. 



TABLE IV: Symmetry classes that support non trivial charge 
or fermion parity pumping cycles. 



an integer number of charges is transported towards or 
away from the point defect. Classes BDI and D define 
a fermion parity pump. We will discuss these two cases 
separately. 

We note in passing that the J = column of Table U 
also applies to topological textures^ for which d — D. For 
example a spatially dependent three dimensional band 
structure 'H(k, r) can have topological textures analogous 
to Skyrmions in a 2D magnet. Such textures have re- 
cently been analyzed by Ran, Hosur and Vishwanath^^ 
for the case of class D, where they showed that the Z2 
invariant characterizing the texture corresponds to the 
fermion parity associated with the texture. Thus, non 
trivial textures are fermions. 



B. Class D, BDI: Fermion Parity Pump 

Adiabatic cycles of point defects in class D and BDI 
are characterized by a Z2 topological invariant. In this 
section we will argue that a non trivial pumping cycle 
transfers a unit of fermion parity to the point defect. 
This is intimately related to the Ising non-Abelian statis- 
tics associated with defects supporting Majorana bound 
states. 

Like the point defect in class Dili {5 = 1), the temporal 
pump (6 — 0) in class D occupies an entry in Table |I] 
similar to the line defect [5 = 2) in class All, so we 
expect a formula that is similar to the formula for the 
2D quantum spin Hall insulator. This is indeed the case, 
though the situation is slightly more complicated. The 
Hamiltonian H(k, r, is defined on a base space T'^ x 
S'^~^ X S^. In appendix IfI we will show that the invariant 
can be written in a form that resembles p.25p . 



A. Class A, AI, All: Thouless Charge Pumps 

The integer topological invariant characterizing a 
pumping cycle in class A is simply the Chern number 
characterizing the Hamiltonian ^-L{\^,Y,t)r^^^. Imposing 
time reversal symmetry has only a minor effect on this. 
For 9^ = — 1 (Class AH), an odd Chern number violates 
time reversal symmetry, so that only even Chern num- 
bers are allowed. This means that the pumping cycle 
can only pump Kramers pairs of electrons. For <d^ = +1 
(Class AI) all Chern numbers are consistent with time 
reversal symmetry. 

The simplest charge pump is the ID model introduced 
by Thouless^. A continuum version of this model can 
be written in the form, 

H(fc,i) = vka^ + {mi{t) + ek'^)a.^ + m2it)ay. (5.1) 

When the masses undergo a cycle such that the phase of 
nil + i'm.2 a single electron is transmitted down the wire. 
In this case, Ji^k^t) has a non zero first Chern number. 
The change in the charge associated with a point in a ID 
system is given by the difference in the Chern numbers 
associated with either side of the point. Thus, after a 
cycle a charge e accumulates at the end of a Thouless 
pump. 

A two dimensional version of the charge pump can be 
developed based on Laughlin's argument!^ for the integer 
quantum Hall effect. Consider a 2D v = 1 integer quan- 
tum Hall state and change the magnetic flux threading a 
hole from to h/e. In the process, a charge e is pumped 



(i!(27r) 



Tr(J-'^ 



Tl/2 



dTi/2 



mod 2, 



(5.2) 

where 7i/2 is half of the base manifold, say fci G [0,7r], 
and the Chern-Simons form Q2d-i is generated by a con- 
tinuous valence frame Wt,(k, r, t)|fei=o,7r that obeys certain 
particle-hole gauge constraint. This is more subtle than 
the time reversal gauge condition p.24p for line defects 
in All and point defects in DHL Unlike p.24|) . we do not 
have a computational way of checking whether or not a 
given frame satisfies the constraint. Nevertheless, it can 
be defined, And in certain simple examples, the particle- 
hole constraint is automatically satisfied. 

The origin of the difficulty is that unlike time reversal 
symmetry, particle hole symmetry connects the conduc- 
tion and valence bands. The gauge constraint therefore 
involves both. Valence and conduction frames can be 
combined to form a unitary matrix, 



G 



k,r,t 



x„(k,r,i) Wc(k,r,t) e U{2n). (5.3) 



The orthogonality of conduction and valence band states 
implies that 



= 0. 



(5.4) 



In general, we call a frame G : dTi/2 U{2n) particle- 
hole trivial if it can continuously be deformed to a con- 
stant while satisfying (j5.4p throughout the deformation. 
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The Chern Simons term in (j5.2p requires a gauge that is 
built from the valence band part of a particle-hole trivial 
frame. 

Though the subtlety of the gauge condition makes a 
general computation of the invariant difficult, it is possi- 
ble to understand the invariant in the context of specific 
models. Consider, a theory based on a point defect in 
the d dimensional version of (j4.6|) . 



9=0 ^ 



^ .4- 



<p (t) 



(5.5) 



Here F and 7 are 2'^ x 2'' Dirac matrices, and we suppose 
that for fixed t, the d dimensional mass vector 0(r, t) has 
a point topological defect at ro (t) . If Ref. [H we argued 
that adiabatic cycles for such point defects are classified 
by a Pontrjagin invariant similar to that discussed in sec- 
tion |IVE31 This may also be understood in terms of the 
rotation of the "orientation" of the defect. Near the de- 
fect, suppose (j){r,t) = 0{t) ■ (r — ro(i)), where 0{t) is 
a time dependent 0{d) rotation. In the course of the 
cycle, the orientation of the topological defect, charac- 
terized by 0{t) goes through a cycle. Since for d > 3, 
7ri(0((i)) = Z2, there are two classes of cycles. As shown 
in Ref. l55l the non trivial cycle, which corresponds to 
a 27r rotation changes the sign of the Majorana fermion 
wavefunction associated with the topological defect. We 
will argue below that this corresponds to a change in the 
local fermion parity in the vicinity of the defect. For 
d = 2, 7ri(0(2)) = Z. However, the change in the sign of 
the Majorana bound state is given by the parity of the 
0(2) winding number. In theories with more bands, it is 
only this parity that is topologically robust. 

In d ~ 1, the single F matrix in the 2 band model does 
not allow for continuous rotations. Consider instead Ki- 
taev's model^ for a ID topological superconductor with 
at time dependent phase, 

'H{k,t) ^ (tcosfc — ^J')'Tz + Ai(i) sinfcr^ + A2(i) sinfcry. 

(5.6) 

In this case it is possible to apply the formula (|5.2p be- 
cause on the boundary (9T, which is fc = or fc = tt 
the Hamiltonian is independent of t, so that the gauge 
condition (|5.4p is automatically satisfied. Moreover, the 
second term in (|5.2p involving the Chern Simons integral 
is equal to zero, so that the invariant is simply the inte- 
gral of T{x,t) over Ti/2- It is straightforward to check 
that this gives v = 1. 

In order to see why this corresponds to a pump for 
fermion parity, suppose a topological superconductor is 
broken in two places, as shown in Fig. [HI At the ends 
where the superconductor is cut there will be Majorana 
bound states. The pair of bound states associated with 
each cut define two quantum states which differ by the 
parity of the number of electrons. If the two ends are 
weakly coupled by electron tunneling then the pair of 
states will split. Now consider advancing the phase of the 
central superconductor by 27r. As shown in Refs. Issllssl . 
the states interchange as depicted in Fig. [S] The level 
crossing that occurs at tt phase difference is protected by 



(a) 



9=0 



^ ^ 'fermion parity^ _ ^ ^ 



(b) 




(P (t) 



FIG. 9: A one dimensional fermion parity pump based on a 
ID topological superconductor, which has Majorana states at 
its ends. When the phase of the central superconductor is ad- 
vanced by 2-n the fermion parity associated with the pairs of 
Majorana states inside each circle changes. Thus fermion par- 
ity has been pumped from one circle to the other, (b) shows 
the evolution of the energy levels associated with a weakly 
coupled pair of Majorana states as a function of phase. The 
level crossing at = tt is protected by the local conservation 
of fermion parity. 



the conservation of fermion parity. Thus, at the end of 
the cycle, one unit of fermion parity has been transmitted 
from one circled region to the other. 

The pumping of fermion parity also applies to adia- 
batic cycles of point defects in higher dimensions, and 
is deeply connected with the Ising non-Abelian statistics 
associated with those defects^. 



VI. CONCLUSION 

In this paper we developed a unified framework for 
classifying topological defects in insulators and super- 
conductors by considering Bloch/BdG Hamiltonians that 
vary adiabatically with spatial (and/or temporal) pa- 
rameters. This led to a generalization of the bulk- 
boundary correspondence, which identifies protected 
gapless fermion excitations with topological invariants 
characterizing the defect. This leads to a number of ad- 
ditional questions to be addressed in future work. 

The generalized bulk-boundary correspondence has the 
flavor of a mathematical index theorem, which relates an 
analytic index that characterizes the zero modes of a sys- 
tem to a topological index. It would be interesting to see 
a more general formulation of this relationSii^^ that ap- 
plies to the classes without chiral symmetry that have 
Z2 invariants, and goes beyond the adiabatic approxima- 
tion we used in this paper. Though the structure of the 
gapless modes associated with defects make it clear that 
such states are robust in the presence of disorder and in- 
teractions, it would be desirable to have a more general 
formulation of the topological invariants characterizing a 
defect that can be applied to interacting and/or disor- 
dered problems. 
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An important lesson we have learned is that topologi- 
cally protected modes can occur in a context somewhat 
more general than simply boundary modes. This expands 
the possibilities for engineering these states in physical 
systems. It is thus an important future direction to ex- 
plore the possibilities for heterostructures that realize 
topologically protected modes. The simplest version of 
this would be to engineer protected chiral fcrmion modes 
using a magnetic topological insulator. The perfect elec- 
trical transport in such states could have far reaching 
implications at both the fundamental and practical level. 
In addition, it is worth considering the expanded possi- 
bilities for realizing Majorana bound states in supercon- 
ductor heterostructures, which could have implications 
for quantum computing. 

Finally, it will be interesting to generalize these topo- 
logical considerations to describe inherently correlated 
states, such as the Laughlin state. Could a fractional 
quantum Hall edge state arise as a topological line defect 
in a 3D system? Understanding the topological invari- 
ants that would characterize such a defect would lead to 
a deeper understanding of topological states of matter. 
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Appendix A: Periodicity in symmetry and dimension 



In this appendix we will establish the relations (|2.5l2.6p 
between the iC-groups in different position-momentum 
dimensions {D, d) and different symmetry classes s. We 
will do so by starting with an arbitrary Hamiltonian in 
Kf{s]D,d) and then explicitly constructing new Hamil- 
tonians in one higher position or momentum dimension, 
which have a symmetry either added or removed. The 
new Hamiltonians will then belong to Kr{s + 1; D , d + 1) 
or Kf{s — 1;D + l,d). The first step is to identify 
the mappings and show they preserve the group struc- 
ture. This defines group homomorphisms relating the K 
groups. The next step is to show they are isomporphisms 
by showing that the maps have an inverse, up to homo- 
topic equivalence. 



1. Hamiltonian mappings 

There are two classes of mappings: those that add sym- 
metries and those that remove symmetries. These need 
to be considered separately. 

We consider first the symmetry removing mappings 
that send a Hamiltonian "He with chiral symmetry to 
a Hamiltonian Unc without chiral symmetry. Suppose 
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FIG. 10: Hamiltonian mappings (|A1|) and (|A3|) are drawn on 
the left and right clocks respectively. Solid (dotted) arrows 
represents addition of one momentum (spatial) dimension. 




FIG. 11: Suspension E(T'' x S^). The top and bottom of 
the cylinder E(T'* x S°) x [-7r/2,7r/2] are identified to two 
points. 



{?^c(k,r),H} 
define 



0, where H is the chiral operator. Then 



-H„c(k,r,6l) = cos6lHc(k,r) + sin6in (Al) 

for — 7r/2 < < tt/2. This has the property that at 
6 — ±7r/2 the new Hamiltonian is ±H, independent of k 
and r. Thus, at each of these points we may consider the 
base space T'^ x defined by k and r to be contracted 
to a point. The new Hamiltonian is then defined on the 
suspension I](T'^ x S^) of the original base space (see 
fig fTTj) . If we treat the original base space as a. d + D 
dimensional sphere, then the suspension is a d + D + 1 
dimensional sphere. 

Without loss of generality we assume He is flattened, 
so that nl = l. Since {He, H} = it follows that Hl^ = 
1 as well. The second term in (|A1[) violates the chiral 
symmetry. Thus, if He belongs to the complex class AIII 
(with no anti unitary symmetries), then Hnc belongs to 
class A. Eq. (|Al[) thus provides a mapping from class AIII 
to class A. 

For the real classes, which have anti unitary sym- 
metries, the second term will violate either particle- 
hole symmetry or time reversal symmetry, depending on 
whether is a momentum or position type variable (odd 
or even under and 2) . This will lead to a new non chi- 
ral symmetry class related to the original class by either 
a clockwise or counter clockwise turn on the symmetry 
clock (figflU]). To determine which it is, note that if we re- 
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quire [e,S] = then (eS)^ = G^S^ = (-l)(«-i)/2. The 
unitary chiral symmetry operator (satisfying = 1) can 
then be written 

n = f'-^^/'^QE. (A2) 

It foUows that if 9 is momentum like, then time reversal 
symmetry is violated when s = 1 mod 4, while particle 
hole is violated when s = 3 mod 4. This corresponds 
to corresponds to a clockwise rotation on the symmetry 
clock, s s + 1. If 6' is position like then s ^ s — 1. 

We next build a chiral Hamiltonian from a non chi- 
ral one by adding a symmetry. This is accomplished by 
doubling the number of bands in a manner similar to the 
doubling employed in the Bogoliubov de Gennes descrip- 
tion of a superconductor. We thus write 

'Hc{k,r,e) = cos6''H„c(k,r) (8)T2 +smei(^Ta, (A3) 

where a = x ov y. Here r are Pauli matrices that act on 
the doubled degree of freedom. As in (|Aip . (jA3p gives a 
new Hamiltonian defined on a base space that is the sus- 
pension of the original base space. If H^c = 1 it follows 
that nl = I, so the energy gap is preserved. It is also 
clear that the new Hamiltonian has a chiral symmetry 
because it anticommutcs with H = iTzTq. Thus, liT-Lnc is 
in class A, then "He is in class AIII. 

For the real symmetry classes a = x ov y must be cho- 
sen so that the second term in (jA3[) preserves the original 
anti unitary symmetry of "Hnc- This depends on the orig- 
inal anti unitary symmetry and whether 6 is chosen to be 
a momentum or a position variable. For example, if T-Lnc 
has time reversal symmetry, O, and is a momentum 
(position) variable, then we require a = y [a = x). In 
this case, T-Lc has the additional particle-hole symmetry 
S = T^Q (S = iTyQ) that satisfies = 9^ (S^ = -G^). 
A similar analysis when "Hnc has particle hole symme- 
try allows us to conclude that the symmetry class of He 
is given by a clockwise rotation on the symmetry clock, 
s — > s -|- 1, when is a momentum variable. When 9 is 
a position variable, s — > s — 1 gives a counter clockwise 
rotation. 

Equations (|A1[) and (jA3p map a Hamiltonian into a 
new Hamiltonian in a different dimension and different 
symmetry class. It is clear that two Hamiltonians that 
are topologically equivalent will be mapped to topologi- 
cally equivalent Hamiltonians, since the mapping can be 
done continuously on a smooth interpolation between the 
original Hamiltonians. Thus, (jAllA3|) define a mapping 
between equivalence classes of Hamiltonians. Moreover, 
since the direct sum of two Hamiltonians is mapped to 
the direct sum of the new Hamiltonians, the group prop- 
erty of the equivalence classes is preserved. (jAllA3|) thus 
define a /sT-group homomorphism, 

Kr{s;D,d) Kw{s + I; D,d + 1), (A4) 
Kr{s;D,d) —> Kr{s ~ I; D + l,d), (A5) 



2. Invertibility 

In order to establish that (|A4IA5p are isomorphisms 
we need to show that there exists an inverse. This is 
not true of the Hamiltonian mappings. A general Hamil- 
tonian cannot be built from a lower dimensional Hamil- 
tonian using (jAllA3|) . However, we will argue that it 
is possible to continuously deform any Hamiltonian into 
the form given by (|A1[) or (jA3p . Thus, the mappings be- 
tween equivalence classes have an inverse. To show this 
we will use a mathematical method borrowed from Morse 
theory^. 

Without loss of generality we again consider flattened 
Hamiltonians having equal number of conduction and va- 
lence bands with energies ±1. Consider H(k, r, 9), where 
9 S [— 7r/2,7r/2] is either a position or momentum vari- 
able and H is independent of k and r a.t 9 — ±7r/2. 
We wish to show that H(k, r, 9) can be continuously de- 
formed into the form (|A1[) or (|A3[) . To do so we define 
an artificial "action" 

S[Hik,r,9)]= J d9d'^kd^rTT{de'Hdgn). (A6) 

S can be interpreted as a "height" function in the space 
of gapped symmetry preserving Hamiltonians. Given any 
Hamiltonian there is always a downhill direction. These 
downhill vectors can then be integrated into a deforma- 
tion trajectory. Since the action is positive definite, it is 
bounded below. The deformation trajectory must end at 
a Hamiltonian that locally minimizes the action. 

Under the flatness constraint = 1, minimal Hamil- 
tonians satisfy the Euler-Lagrange equation 

d^H + H^0. (A7) 

The solutions must be a linear combination of sin 9 and 
cos 9. The coefficient of sin6' must be constant because 
the base space is compactified to points at6' = ±7r/2. A 
minimal Hamiltonian thus has the form 

H(k,r,6l) =cos6'Hi(k,r) +sine'-Ho. (A8) 

The constraint H(k,r,9)^ = 1 requires 

H^=Hi(k,r)2 = l, {Ho,1^i(k,r)}==0. (A9) 

If ?^(k, r, 9) is non-chiral, then ea. (|A8l) is already in the 
form of cq.dM]) with H = 7^0 and 'Hc(k,r) 'Hi(k,r). 
Til automatically has chiral symmetry due to Eq. (|A9p . 
This shows that (jA4p and (jASP are invertible when s is 
odd. 

If H(k,r,9) is chiral, then both Hq and ■Hi(k, r) an- 
ticommute with the chiral symmetry operator H. Re- 
name Hq =^ Ta and H = iTzTa, where a — x (a ~ y) 
when is a position (momentum) variable. It follows 
that {HijTx} — {HijTy} = 0, so we can write 



for F = M, C. 



'Hi(k,r) = h{k,r) r^. 



(AlO) 
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Eq. 



thus takes the form of Ea. (|A3l) with % 



Since anti-commutes with either or S, /i(k, r) carries 
exactly one anti-unitary symmetry and is therefore non- 
chiraL This shows that (IA4I) and (IA5I) are invertible when 



s IS even. 



Appendix B: Representative Hamiltonians, and 
classification by winding numbers 

In this appendix we construct representative Hamilto- 
nians for each of the symmetry classes that are built as 
linear combinations of Clifford algebra generators that 
can be represented as anticommuting Dirac matrices. 
This allows us to relate the integer topological invariants, 
corresponding to the Z and 2Z entries in Table I, to the 
winding degree in maps between spheres. Similar con- 
struction for defectless bulk Hamiltonians can be found 
in Ref.frU by Ryu, et.al. In general, Hamiltonians do 
not have this specific form. However, since each topolog- 
ical class of Hamiltonians includes representatives of this 
form, it is always possible to smoothly deform 'H(k, r) 
into this form. 

The simplest example of this approach is the familiar 
case of a two dimensional Hamiltonian with no symme- 
tries (class A). A topologically non trivial Hamiltonian 
can be represented as a 2 x 2 matrix that can be ex- 
pressed in terms of Pauli matrices as H(k) = h(k) • a. 
The Hamiltonian can then be associated with a unit vec- 
tor d(k) = h(k)/|h(k)| G S'^. It is then well known that 
the Chern number characterizing 'H(k) in two dimensions 
is related to the degree, or winding number, of the map- 
ping from k to 5^. This approach also applies to higher 
Chern numbers characterizing Hamiltonians in even di- 
mensions d = 2n. In this Hamiltonian that is a 
combination of 2ri, -|- 1 2" x 2" Dirac matrices, and can 
be associated with a unit vector d € 5^" . 

For the complex chiral class AIII, the U{n) winding 
number characterizing a family of Hamiltonians can sim- 
ilarly be expressed as a winding number on spheres. For 
example, in d = 1, a chiral Hamiltonian can be written 
'H(fc) = hx{k)ax + hy{k)ay (so {H, a^} = 0), and is char- 
acterized by d(fc) € S^. The integer topological invariant 
can then be expressed by the winding number of d(fc). 
Similar considerations apply to the integer invariants for 
chiral Hamiltonians in higher odd dimensions. 

For the real symmmetry classes we introduce "position 
type" Dirac matrices and "momentum type" Dirac 
matrices 7^. These satisfy {F^,Fy} = 26^^, {7^,7^} = 
25ij, {r^,7j} = and are distinguished by their sym- 
mmetry under anti unitary symmetries. If there is time 
reversal symmetry we require 



[F^,e] = {7.,e} = o, 

while with particle-hole symmetry, 

{F,„S} = [7„S]=0. 



(Bl) 



(B2) 
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TABLE V: Examples of Dirac matrices for (p,q) = (s,0). 



For a Hamiltonian that is a combination of p momen- 
tum like matrices 71,..., p and <? + 1 position like matrices 
To,...,?, 



H(k,r) =R(k,r)-r + K(k,r)-7, 

the coefficients must satisfy the involution 

R(-k,r) = R(k,r), 
K(-k,r) = -K(k,r). 

This can be characterized by a unit vector 



d(k,r) 



(K,R) 



IRI- 



(B3) 



(B4) 
(B5) 



(B6) 



where S^'^'^ is a (p-l-q)-sphere in which p of the dimensions 
are odd under the involution (jB5|) . 

The symmetry class s of 'H(k, r) is related to the in- 
dices ij), q) characterizing the numbers of Dirac matrices 
by 



mod 



(B7) 



To see this, start with a Hamiltonian "Hq — -Ro(k, r)ro 
that involves a single 1x1 position like "Dirac matrix" 
Fq = 1, so (p, (?) = (0,0). This clearly has time rever- 
sal symmetry, with Q = K , and corresponds to class AI 
with s = 0. Next, generate Hamiltonians "Hs with dif- 
ferent symmetries s by using the Hamiltonian mappings 
introduced in appendix IA II Both the mappings (jAl[) and 
(|A3|) define a new Clifford algebra with one extra genera- 
tor that is either position or momentum type. The map- 
pings that correspond to clockwise rotations on the sym- 
metry clock (s — > s + 1) introduce an additional position 
like generator (p — > p + 1), while the mappings that cor- 
respond to counterclockwise rotations (s — > s — 1) intro- 
duce an additional momentum like generator (g — g-f-l). 
Eq. ()B7P follows because this procedure can be repeated 
to generate Hamiltonians with any indices {p,q). Some 
examples are listed in table |V] 

The integer topological invariants in Table U (which 
occur when s — 5 is even) can be related to the winding 
degree of the maps d : S^^'^ — ?> 5^'+*. This can be non 
zero when the spheres have the same total dimensions. In 
light of (|B7|) . (p, q) can always be chosen so that d + D = 
p + q. The anti unitary symmetries impose constraints 
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on the possible values of these winding numbers, which 
depend on the relation between 5 = d-- D and s = p — q. 

The involutions on S'^'^^ and 5'^'*'' have opposite ori- 
entations when (5 — s = 2 or 6 mod 8, and therefore an in- 
volution preserving map S'^'^^ — > 5*^+' can have non-zero 
winding degree only when S — s = or 4 mod 8. Symme- 
try gives a further constraint on the latter case. Consider 
a sphere map Sg ^ — ^, where the involutions on the 
spheres send (6*, </>) i~7> (6*, </> -I- tt) and '-^ {''^t'^)- In 

order for f{9, 0) = f{0, </)-!- tt), the winding number must 
be even. Together, these show 



deg e 




for (5 — s = (mod 
for S — s = A (mod 
otherwise 



(B8) 



This gives a topological understanding of the Z's and 
2Z's on the periodic table in terms of winding number, 
which can be identified with the more general analytic 
invariants, namely Chern numbers for non-chiral classes 



1 



Tdy,SD 



Tr LF^ 



(B9) 



and winding numbers of the chiral flipping operator 
(7(k, r) for chiral ones (See Eqs. 14.11 and 142)) ). 



{d + D)\{2m)- 



Tr((qrfgt) 



r\d+D\ 



(BIO) 



The Z2 's on the periodic table are not directly charac- 
terized by winding degree, but rather through dimen- 
sional reduction. Given a Hamiltonian H(k, fci, fc2, r) 
with s = S mod 8, its winding degree mod 2 de- 
termines the Z2-classification of its equatorial offspring 
'Hfc2=o(k, fci, r) and "Hfci 2=o(kj r). For example, topolog- 
ical insulators in two and three dimensions are equatorial 
restrictions of a four dimensional model d : 5*^ ^ 5*^ 
with unit winding number. Around the north pole, the 
Hamiltonian has the form 

n{k, k^) = {m + ek'^)ni + k • /iga + ki^i2 (Bll) 

and on the equator fc4 = 0, this gives a three dimen- 
sional Dirac theory of mass m around k = that locally 
describes 3D topological insulators Bi2Sc3 and Bi2Te3 
around F. 



Appendix C: Zero Modes in the Harmonic Oscillator 
Model 

We present exact solvable soliton states of Dirac-type 
defect Hamiltonians. These include zero modes at a point 
defect of a Hamiltonian in the chiral class AHI, and chiral 
modes along a line defect of a Hamiltonian in the non- 
chiral class A. We establish the connection between the 
two kinds of boundary modes through the Hamiltonian 
mapping (|A1I) . 




FIG. 12: Energy spectrum of Hamiltonian (|C4|) . The zero 
mode of positive chirality at fcy = corresponds the chi- 
ral mode that generates the mid-gap /cy -linear energy spec- 
trum. 



A non-trivial chiral Hamiltonian isotropic around a 
point defect at r = is a Dirac operator 



-17 • V -I- r • F 



(CI) 



where the chiral operator is H = 11^=175^^1 ^nd its 
adiabatic limit e~^^''^He^^''^ = k-7 + r-r has unit winding 
degree on S'^d-i {(k,r) : fc^ + = 1}. 



(C2) 



and the spectrum is determined by the quantum numbers 
Uj > of the harmonic oscillator and the parities S,j of 
the mutually commuting matrices ijjTj, for j ~ I, . . . ,d. 



(C3) 



The unique zero energy state |5'o), indexed by Uj = 
and =1, has positive chirality H = -1-1, and is expo- 
nentially localized at the point defect as 4*0 (r) (x e^^^ . 

Next we consider a non-chiral Hamiltonian isotropic 
along a line defect. 



n{kn) = fc||H-i7-V + r-F 



(C4) 



where fcy is parallel to the defect line, r and V are 
normal position and derivative. Its adiabatic limit 
e-''^"n{k^^)e'^-'' = /s||n + k • 7 + r • f is related to that 
of ()C1|) by (l,l)-pcriodicity, and has unit winding degree 
on S^'^ = {(fc||,k,r) : k^^'+ \k\^ + \r\^ = 1}. The zero 
mode l^'o) of (|C1|) gives rise to a positive chiral mode, 

H(fc||)|*o) = fc||n|*o> = +fc|||*o). 

The two examples verified bulk-boundary correspon- 
dence through identifying analytic information of the 
defect-bound solitons and the topology of slowly spatial 
modulated theories far away from the defect. The single 
zero mode of (jCip and spectral flow of (jC4[) are equated 
to unit winding degree of an adiabatic limit. In general, 
bulk-boundary correspondence is mathematically sum- 
marized by index theorems that associate certain analytic 
and topological indices of HamiltonianSi '^^'^^i'^^'^^i^^'^^ 
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Appendix D: Invariant for Point Defects in Class D 
and BDI 



We follow the derivation given in Ref.[55|, which was 
based on Qi, Hughes and Zhang's formulation of the 
topological invariant characterizing a three dimensional 
topological insulator—. For a point defect in d dimen- 
sions, the Hamiltonian "^(14, r) depends on d momentum 
variables and d — 1 position variables. We introduce 
a one parameter deformation ■H(A,k, r) that connects 
H(k,r) at A = to a constant Hamiltonian at A = 1, 
while breaking particle-hole symmetry. The particle-hole 
symmetry can be restored by including a mirror image 
-H(A,k,r) = -S-H(-A,k,r)S-i for -1 < A < 0. For 
A = ±, (k, r) can be replaced by a single point, so the 2d 
parameter space (A, k, r) is the suspension E(r'* x S"^^^) 
of the original space. The Hamiltonian defined on this 
space is characterized by its d'th Chern character 




FIG. 13: (a) Schematic of the base space T^''-^ = T'' x S"'-^ 
split into two halves, (b) Division of 'Y^'^~^ into patches A 
and B, each closed under TR and has an individual valence 
frame \um^) that satisfies the gauge condition (|E2|) . 



Appendix E: Invariant for line defects in class All 
and point defects in Dili 



1 

d! V 27r 



S(T''xS''-i) 



Tr[J-'' 



(Dl) 



Due to particle-hole symmetry, the contributions from 
the two hemispheres A > 0, A < are equal. Using 
the fact that the integrand is the derivative of the Chern 
Simons form, Tr[J-'^] = dQ2d-i: we can therefore write 



"=d! 



27r 



l2d~l 



(D2) 



As was the case in Refs.Hili, V can be different for 
different deformations 'H(A,k, r). However, particle-hole 
symmetry requires the difference is an even integer. 
Thus, the parity of (jD2| defines the I2 invariant. 

The Chern Simons form Q2d-i can be expressed in 
terms of the connection A via the general formula 



i2d- 



-i=d 



f dtTr[A{tdA + t'^AY~^] (D3) 
Jo 



In the addition of time reversal symmetry 0^ = 1 or 
equivalently a chiral symmetry H = 0H = Tz, a valence 
frame of the BDI Hamiltonian (|4.ip can be chosen to be 



u(k, r) 



J_fq{k,r) 

V2\ -1 



(D4) 



where q is unitary and 1 is the identity matrix. This 
corresponds the Berry connection A = u^du = ^q^dq 
and Chern-Simons form 



!2d-l 



dt 



t f t 



2 V2 



l]C\.[{qUqr-'] 



(D5) 



This equates the winding number of q (|4.2p to the Chern 
Simons invariant (|4.9p . 



We formulate a topological invariant that characterizes 
line defects in class AH in all dimensions that is analo- 
gous to the integral formula invariant characterizing the 
quantum spin Hall insulator introduced in Ref. [tsI . This 
can be applied to weak topological insulators in three 
dimensions with dislocation around a line defect. The 
invariant can be indirectly applied to strong topological 
insulators through decomposition into strong and weak 
components. As a consequence of the Hamiltonian map- 
ping (|Aip that identifies {s = — 2) and (s = 3, 5 = 1), 
this gives a new topological invariant that classified point 
defects in class Dili in all dimensions. 



Line defects in class All 



The base space manifold is T^'* ^ = T"^ x S"^ ^, where 



is the Brillouin zone and S' 



d-2 



is a cylindrical 



neighborhood that wraps around the line defect in real 
space. Divide the base space into two pieces, 7^/2~^ ^^'^ 
its time reversal counterpart (see fig fTHT a)). We will show 
the Z2-invariant 



(d- l)!(27r 



ld-l 



Tr(J-'^-i) 



r- 



1/2 



Q 



2d-3 



9r, 



1/2 



(El) 

topologically classifies line defects in AH, where the 
Chern-Simons form, defined by (|D3p . is generated by the 
Berry connection Amn = (wm(k, r)|dM„(k, r)), and the 
valence frame Um (k, r) satisfies the gauge condition 



.(k,r) 



i(k, r)|0u„(— k, r)) = constant (E2) 



on the boundary (k, r) e 9T^/2 ^ ■ 

The non-triviality of the Z2-invariant is a topological 
obstruction to choosing a global continuous valence frame 
|um(k, r)) that satisfies the gauge condition (|E2p on the 
whole base space 7~2ti-2^ j£ there is no topological ob- 
struction from the bulk^S, the gauge condition forces the 
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valence frame to be singular at two points, depicted in 
figfTHTb). related to each other by time reversal. One re- 
moves the singularity by picking another valence frame 
locally defined on two small balls enclosing the two singu- 
lar points, denoted by B in fig flBl b). We therefore have 

two valence frames |um^^(k, r)) defined on two patches 
of the base space, A = T'^'^~^\B and B, each obeying the 
gauge condition (|E2[) . 

The wavefunctions on the two patches translate into 
each other through transition function 



C::(k,r) = (^.;^(k,r)K(k,r))GC/(fc) 



(E3) 



on the boundary dB sa S'^''-^ y 52^-3 rp^^ function 
behavior on the two disjoint {2d — 3)-spheres is related 
by time reversal. The topology is characterized by the 
winding of t^^ : S^"^'^ -> U (k) on one of the spheres. 



(d-2)! 



(2d - 3)!(27ri) 

{d ~ iy.{2TTiY- 



Tr 



22d-3 0-2(1-3,1 



S2d-3 



(E4) 
(E5) 



The two integrals can be evaluated separately. Since 
dQ2d-3 = Tr(J''^-i), Stokes' theorem teUs us 



[ Tr(J-' 



Combining these into eq. ljESp identifies the Z2-invariant 
(|E1[) with the winding number of the transition function. 
The curvature term in ()E1|) is gauge invariant. Any 



gauge transformation on the boundary dT^^^ ^ respect- 
ing the gauge condition (|E2[) has even winding number 
and would alter the Chern-Simons integral by an even in- 
teger. The gauge condition is therefore essential to make 
the formula non- vacuous. 



a. Spin Chern number 

A quantum spin Hall insulator is characterized by its 
spin Chern number = {n^ — n^)/2. We generalize this 
to time reversal invariant line defects of all dimensions 
by equating it with Eq. (|El[) . This applies in particular 
to a model we considered for a linear Josephson junction 
in section UlIDI 

A spin operator 5 is a unitary operator, square to 
unity, commutes with the Hamiltonian, and anticoni- 
mutes with the time reversal operator. The valence spin 
frame 



K4(k,r)) =e|^.i-„(~k,r)) 



(E6) 



automatically satisfies the time reversal gauge constraint 
(IE2[) . It is straightforward to check that the curvature 




FIG. 14: (a) Schematic of the suspension ET^'^'V (b) De- 
composition into patches A and B. 



and Chern-Simons form can be split as direct sums ac- 
cording to spins. 



J"(k,r) = T^{k,r)®J-^{-k,r) 
4_3(k,r)eQ|,_3( 



(E7) 

Q2rf_3(k,r) = Q|,_3(k,r)eQ|,_3(-k,r)* (E8) 



Again assuming that there is no lower dimensional 
"weak" topology, the f-frame can be defined everywhere 
on 7~2<i-2 -^^ith a singularity at one point, say in 7^^^2~^' 
and the ^-frame is singular only at the time reversal of 
that point. 

The curvature term of (|Eip splits into two terms 



/ Tri^^-')^ I -[ 

'wo L 1/2 



/7-2d-2\7-2d-2 



Tr(^l'-i) 
(E9) 

And the two spin components of the Chern-Simons term 

f-,-7-2d-2 Q2d-3 add up into 
" ' 1/2 



Qld-3 



dTf. 



7-2d-2\7-2< 



Tr(-F' 



(ElO) 



by Stokes theorem. 

Combining these two, we equate (|E1[) to the spin Chern 
number 



■d-l 



Tr(J-.^-i) 



(Ell) 



Time reversal requires ntot = + = and therefore 
n„ = (n^ - n;)/2 = n^. 

2. Point defects in class Dili 

The base space manifold is T^'^^^ =T'^ x S'^'^. The 
Hamiltonian mapping (|A1[) relates a point defect Hamil- 
tonian 'H(k, r) in class Dili to a line defect Hamiltonian 
■H(k, r,0) = cos6''H(k, r) -f sinflH in class AH, where 
H = i9E! is the chiral operator, (k, r, 0) S ST^'*"^ (see 
fig fTW a)') and 9 is odd under time reversal. The line de- 
fect Hamiltonian 'H(k, r, 6) is topologically characterized 
by the generalization of (|Eip . which was proven to be 
identical to the winding number (jE4[) of the transition 
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function t^^ (see fig lMT bl for the definition of patches 
A and B). We wih utihze this to construct a topofogical 
invariant that characterizes point defects in class DHL 

Set O = TyK and 2 = t^K under an appropriate 
choice of basis. A canonical valence frame of 'H(k, r,0) 
can be chosen to be 



(k,r,( 



sm 



f)9(k,r) 
1 



(E12) 



where ^(k, r) £ U (fc) is from the canonical form of the 
chiral Hamiltonian H(k, r) in eq. (|4.ip . 1 is the kxk iden- 
tity matrix, and the valence frame is non-singular every- 
where except at 9 = —tt/2. There is a gauge transfor- 



mation u5 — > w'^ = u^t^^ everywhere except 9 = ±7r/2 



such that the new frame satisfies the gauge condition 
(|E2P i^ A valence frame on patch B can be constructed 
by requiring u^{9) = Qu^{—9) around 9 = —7t/2. 

The Z2-topology is characterized by the evenness or 
oddness of the winding number of t^^ as in (jE4[) . This 
can be evaluated by the integral along the equator = 



(d-iy. 



{2d ~ iy.{2TTi)d' JTdy,s'i-i 



Tr 



2d-l 



where = u^t^^ is a solution to the gauge condition 



(E13) 
idition 

(E14) 



where the constant in eq.(|E2j) is chosen to be icjy. 

The winding number (|E13I) can also be expressed as a 
Chern-Simons integral. 



E2p . or equivalently t^^ satisfies 

g(k,r)=t^^(-k,rKi^^(k,r)^ 



d!(27r) 



(E15) 



where Q^^^ arc the Chern-Simons form generated by 
valence frames u^l^ . Restricted to 6* = 0, ()E12|) gives 
w^(k,r) = ^((7(k,r),-l). Following dHl]), the first 

term of (|E15[) equals half of the winding number of 
which is guaranteed to be zero by time reversal and 
particle-hole symmetries. And therefore point defects in 
Dili are classified by the Chern-Simons invariant 



V = — . 
d\ \ 2-K 



Q,2d-\ mod 2 



(E16) 



where the Chern-Simons form is generated by a valence 
frame that satisfies the time reversal gauge constraint 
(|E2l). 

Note that the integrality of the Chern-Simons integral 
(jE16p is a result of particle-hole symmetry. Forgetting 
time reversal symmetry, point defects in class D are clas- 
sified by the Chern-Simons invariant v ~ 2v (jD2|) with 
a factor of 2. Time reversal symmetry requires the zero 
modes to form Kramers doublets, and therefore v ^ 2v 
must be even. A gauge transformation in general can al- 
ter V by any integer. Thus, similar to the formula in class 
All, the time reversal gauge constraint (|E2p is essential 
so that (jE16[) is non- vacuous. 





] 




-^1/2 












\ <!>/ 
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■^1/2 



deformation 



FIG. 15: Deformation of Hamiltonian and its base manifold 
so that the boundary dTi/2 is shrink to a point. 



a. Fixed points formula in ID 



Tiiqkdql) (E17) 



We here identify (|E16[) . or equivalently (IE13I). to a 
fixed point invariant in 1 dimension. In ref.[82|, Qi, 
Hughes and Zhang showed that ID TRI superconductors 
arc Z2-classificd by the topological invariant 

(-lf = ^fe4expfi 
^ ^ Pf(<Zfc=o) V2 

under the basis Q = TyK and S ~ t^K, where qk is 
the chiral flipping operator in (|4.1|) . Time reversal and 
particle- hole symmetry requires qk = ^Q-k- Hence the 
Pfaffians are well defined as qk is antisymmetric at the 
fixed points = 0,7r. 

Using the gauge condition ()E14[I , we can expressed the 
Pfaffians as Pi{Qqk=o,Tr) = det(ife=o,7r)Pf(o'y), where t^^ 
is abbreviated to tk- 



exp 



Tr ( tkdtl 



(E18) 



Pf(9fe=o) V Jo 

Substitute (|E14p into the Cartan form Ti (qkdql,) gives 

(E19) 



Tr qkdq 



Tr ( t^kdt 



Tr tkdt 



Combining these into (jE17p . 
(-1)' 



exp 



exp 



1 

2 JO 
1 
2 



Tr U 



Tr 



fcrfil,-M4))(E20) 
tkdtl)) (E21) 



which agrees (|E13|) . 



Appendix F: Invariant for fermion parity pumps 



In the appendix, we will show that a Z2-invariant ()5.2|) 
under a particlc-holc gauge constraint topologically clas- 
sified fermion parity pumps in dimension (5 = and class 
D or BDI. (See sec lVBj for the fuh statement.) 

We will show (|5.2p using a construction similar to a 
reasoning in Moore and Balents Q . We consider a defor- 
mation of the Hamiltonian along with the base manifold, 
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so that the boundary 97i/2 is deformed into a single point 
(see fig [T5|) . Let s G [0,1] be the deformation variable, 
and denote 7^^2('') ^^'^ dTi/2{s) be the corresponding 
deformation slices. Chern invariant 



Tr(J-'^) 



(Fl) 



integrally classifies Hamiltonians on the half manifold 
^1/2^^ ~ Pa-rticle-hole symmetry requires opposite 
Chern invariant on the other half. A different choice of 
deformation could only change the Chern invariant by 
an even integeri^ Hence the Chern invariant modulo 2 
defines a Z2-invariant. 



The Chern integral can be further deformed and de- 
composed into 

/ Tr(J-'^)+ f ds I Tr(J-'*) 

Jt+{s=o) Jo JdT+ [s) 



I Trl.F'^) - / 



Q 



2d-l 



ar+2(s=o) 



(F2) 



where Stokes' theorem is used and the negative sign is 
from a change of orientation of the boundary. This proves 
(|5.2p . The particle- hole gauge constraint is built-in since 
the Gk.r,t(s) deforms into a constant at s = 1 while re- 
specting particle-hole symmetry (|5.4|) at all s. 
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Bulk topology is a lower dimenional topological obstruc- 
tion coming purely from the momentum part T'' of the base 
space T"^ x S"*"^. One could mathematically remove the 
bulk topological obstruction by adding a defectless Hamil- 
tonian •^(k, r) = H(k, r) Ho(k). 

Again, assume the bulk has trivial topology. Or other- 
wise remove the topological obstruction mathematically by 
adding a defectless Hamiltonian. 

Deformation of Hamiltonians on 97i/2(s), for s G [0, 1] has 
dimension 5 — (d— 1) — (d— 1 + 1 + 1) = —2, and is classified 
by 2Z. 



